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The  development  of  the  gasdynamic  laser  is  one  of  the  more  important 

developments  of  the  past  decade.  This  laser,  based  on  early  works  of  Patel 

(1  2) 

and  others,  ’  obtains  population  inversion  leading  to  laser  action  in  C02  -  N2 
mixtures  by  rapid  expansion  of  the  gases  to  supersonic  speeds.  Relatively  high 
powers,  0(60Kw),  in  continuous  wave  operation  have  been  achieved  by  this  method, 
however,  at  very  limited  efficiency,  G(h  ~  IS) ,  based  on  radiated  energy  to  flow- 
enthalpy.  It  is  the  purpose  of  the  current  program  to  investigate  methods  for 

increasing  this  level  of  efficiency,  the  ideas  of  which  may  find  application  to 
other  laser  systems. 


The  main  objective  of  the  program  is  to  utilize  a  larger  fraction  of  the 
total  enthalpy  of  the  laser  fluid  for  lasing  action  in  an  open  cycle  system. 

The  most  straight-forward  way  to  accomplish  this  is  to  efficiently  expand  the 
gas,  lase,  recompress  and  repeat  this  process  as  many  times  as  possible.  The 
difficulties  to  this  approach  limit  the  viability  of  this  idea  and  form  the 
basis  of  our  research;  they  are:  (1)  low  radiated  energy /flow  energy 
which  makes  the  cycle  very  sensitive  to  component  inefficiency, 

*:2)  largr‘  entr°Py  Production  during  lasing  in  a  channel,  resulting  in  large  loss 
of  total  pressure,  and  (3)  difficulties  during  recompression  (diffuser)  stage 

with  Viscous  effects  that  prevent  efficient  recompression. 

/. 


Increasing  the  radiated  energy /mass  of  fluid  would  solve  problem  area  (1). 
Ope  rat  on  of  a  gasdynamic.  laser  at  higher  temperatures  would  give  this  increased 
effect  and  accordingly  during  this  contract  period,  our  program  investigated  the 
basic  kinetic  processes  of  the  laser  raixt^e  under  high  temperature  conditions. 
Gain  measurements  of  gasdynamic  laser  mixtures  have  been  carried  out  at  stagna¬ 
tion  temperatures  from  1700°K  to  3700°K.  Gain  has  been  observed  over  this  entire 
temperature  range  for  low  C02  concentration  laser  mixtures.  Surprisingly  though, 
our  measurements  show  significantly  higher  gains  than  have  been  predicted  over  a 
wide  range  of  conditions:  It  is  believed  that  this  apparent  gain  enhancement  is 
possibly  due  to  near  resonant  spectral  line  overlap  of  highly  excited  states. 

The  results  of  the  experimental  gain  study  indicates  that  C02  -  N2  mixtures  can 
be  operated  under  high  enthalpy  conditions  if  a  practical  method  of  gas  utiliza¬ 
tion  at  high  temperatures  can  be  found. 

In  addition  we  have  just  begun  a  numerical  investigation  of  the  lasing 
process  using  computer  simulation.  The  primary  aim  is  to  thoroughly  investigate 
the  possibilities  of  reducing  entropy  production  by  lasing  in  an  expanding  flow 
rather  than  a  constant  area  duct.  The  full  system  kinetics  are  being  modeled 
(a  four  temperature  and  vibrational  mode  approach).  One  dimensional  flow  is  as¬ 
sumed  for  simplicity,  and  an  equivalent  one  dimensional  radiation  extraction 
process  based  on  a  laser  amplifier  approach  is  the  preliminary  program.  Initial 
nrn-optimizer  results  tend  to  confirm  early  analytic  predictions  based  on  a 
model  which  ignored  the  complete  kinetic  processes. 

As  pointed  out  in  the  proposal,  a  completely  supersonic  laser  cycle 
(nozzles  as  well  as  diffusers)  may  be  a  more  efficient  device  over  the  more 
conventional  de  Laval  nozzle  (subsonic  and  supersonic,  i.e.  a  mixed  flow)  system. 

It  is  hoped  that  avoiding  transonic  and  subsonic  flow  in  the  diffuser  will  emit 


a  more  efficient  design  than  is  currently  possible  while  scill  carrying  out 
the  function  of  re-energizing  the  vibrational  levels  of  the  tystem.  Part  of  the 
program  is  thus  endeavoring  to  understand  the  non-ideal  aerodynamics  and  kinetics 
of  the  all  supersonic  flow  approach  using  the  Prandtl  Meyer  turning  vanes.  Gain 
measurements  in  the  Prandtl  Meyer  expansion  flow  have  been  completed.  The 
measured  gains  were  always  somewhat  lower  than  predicted,  bit  this  has  been 
traced  to  insufficient  testing  time  in  the  three-inch  shock  tube  facility  so 
that  quasi-steady  flow  in  the  cavity  region  was  not  fully  achieved.  A  boundary 
layer  computer  program  incorporating  the  possibility  of  external  pressure  gradi¬ 
ents  and  non-adiabatic  walls  is  now  beginning  to  give  results.  Preliminary  out¬ 
puts  are  now  available  to  study  their  effect  on  diffuser  performance. 

These  studies  are  individually  explored  in  more  detail  in  the  following 
sections  of  the  report. 

II.  HIGH  ENTHALPY  GAIN  MEASUREMENTS 

Most  gas dynamic  laser  systems  typically  extract  only  -  1%  of  the  total 
flow  energy  (enthalpy)  as  coherent  radiation.  One  of  the  possible  reasons  for 
such  low  performance  levels  in  GDL's  is  that  low  stagnation  temperatures  of  from 
1200°K  to  1600°K  are  commonly  used.  At  such  low  temperatures  the  vibrational 
energy  levels  of  molecular  nitrogen  are  not  fully  excited.  Figure  1  illustrates 
the  amount  of  energy /mole  stored  in  the  vibrational  energy  levels  of  N2  as  a 
function  of  temperature.  The  total  vibrational  energy  increases  steadily  with 
temperature.  While  the  vibratio.ial  energy  in  the  first  level  reaches  a  maximum 
at  about  5000%  the  total  vibrational  energy  of  the  system  greatly  exceeds  the 
energy  stored  in  the  first  excited  vibrational  level  above  2000°.  This  first 
vibrational  level  is  the  level  that  pumps  the  C02  molecule  to  its  upper  laser 
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level.  The  possibility  of  the  upper  vibrational  levels  cascading  down  to  the 
first  excited  vibrational  level  of  N2  in  a  rapidly  expanding  nozzle  flow  „c  Id 
greatly  enhance  the  pumping  action  of  N2  on  the  upper  laser  level  of  C02  and 
hence  convert  more  of  the  flow  energy  of  the  gas  into  coherent  radiation; 

The  use  of  high  temperatures  would  increase  the  relaxation  rate  processes 
in  rapidly  expanding  nozzle  flow  and  hence  hinder  the  freezing  process  that  is 
necessary  to  obtain  a  population  inversion  in  the  nitrogen  and  in  the  upper 
laser  level  of  CO,,.  Although  nitrogen  by  itself  is  easily  frozen,  the  amount 
of  C02  in  the  mixture  has  a  large  effect  on  this  freezing  j recess  and  can  ad¬ 
versely  affect  the  attainment  of  a  population  inversion.  However,  through  the 
use  of  starved  CO.,  mixtures,  the  adverse  effects  of  CO.,  on  the  freezing  process 
can  be  significantly  red„ced«>  Such  starved  mixtures  have  proved  very  effective 
in  GDI's  using  very  high  stagnation  pressures  (po  ,  100  atm>!4>  Also  at  high 
temperatures,  a  CO,,  -  llj  gas  mixture  may  effectively  starve  itself  by  dissoci¬ 
ating  the  C02.  It  is  expected  that  the  products  of  such  a  chemically  dissoci¬ 
ating  flow  (primarily  CO,  0,  0.,,  NO)  would  have  a  smaller  adverse  effect  on  the 
freezing  process  than  CO.,.  Thus,  if  the  CO,  can  be  kept  dissociated  during  the 
earl,  part  of  the  expansion  process,  the  freezing  of  the  11,  vibrational  energy 
states  will  be  greatly  enhanced.  Furthermore,  the  -ossibility  of  C02  being 

formed  m  tl.«  excited  upper  laser  level  upon  recombination  may  add  to  the 
inversion. 

With  the  abo-' e  thoughts  in  mind,  a  series  of  gain  measurements  were  under¬ 
taken  to  uetermine  if  any  of  the  above  effects  of  high  stagnation  temperatures 
would  have  a  favorable  effect  on  the  performance  of  a  GDL.  The  Aerospace 
Research  Laboratory's  7"  I.D.  high  pressure  shock  tube  was  used  to  drive  a  two 
dimensional  sup  rsonic  nozzle.  Gam  measurements  were  made  at  about  the  Mach  5 
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station  of  the  nozzle  using  three  CO^  -  -  He  gas  mixtures  containing  approxi¬ 

mately  10%,  t%,  and  2.5%  C02<  Stagnation  temperatures  were  varied  from  2000°K 
to  5000°K  and  stagnation  pressure  was  kept  constant  at  about  4C  a tins.  The 
experimental  results  indicate  that  high  temperature  effects  may  indeed  have  a 
beneficial  effect  on  the  performance  of  a  GDL.  Detailed  considerations  of  these 
effects  are  discussed  further. 

Experimental  Apparatus  and  Technique 

The  genet al  experimental  technique  used  to  measure  gain  was  essentially 
that  of  Reference  4.  A  shock  tube  is  used  to  produce  a  high  pressure,  high 
temperature  C02  -  N2  -  He  gas  reservoir  behind  a  reflected  shock  wave.  The  gas 
processed  by  the  reflected  shock  is  expanded  through  a  de  Laval  nozzle  producing 
a  population  inversion.  Laser  gain  is  then  measured  at  the  Mach  5  station  of  the 
nozzle.  The  shock  tube  is  adequately  described  in  Reference  4.  The  laser 
optical  system  and  nozzle  design  were  significantly  changed  from  that  of 
Reference  4  in  an  effort  :o  increase  the  signal  to  noise  ratio  and  to  reduce 
damage  to  the  nozzle  blocks  and  germanium  windows.  Acoustic  noise  and  mechanical 
vibrations  and  stress  were  severe  for  this  series  of  tests. 

The  three  gas  mixtures  used  for  these  tests  are  as  follows: 

(1)  10.73%  C02  -  41.02%  N2  -  48.25%  He 

(2)  5.18%  C02  -  47.67%  N2  -  47.15%  He 

(3)  2.63%  C02  -  53.97%  N  -  43.40%  He  . 

These  gas  fixtures  were  obtained  from  Liquid  Carbonic  Corporation  and  are 
pure  grade,  i.e.  99.99 %  pure  or  better. 

A  sketch  op  the  experimental  apparatus  is  shown  in  Figure  2.  The  two 
dimensional  nozzle  blocks  were  machined  from  304  stainless  steel  and  mounted  in 
an  aluminum  flange  that  mated  with  ARL's  7"  I.D.  high  pressure  shock  tube.  This 
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nozzle  test  station  was  designed  to  withstand  impulse  loaded  pressures  of  4000  psi. 
The  nozzle  had  throat  dimensions  <h*)  of  0.1  cm  high  and  15.24  cm  wide,  a  subsonic 
converging  section  with  a  0.4  cm  radius  contour  leading  smoothly  to  the  throat, 
and  supersonic  diverging  section  with  a  half  angle  of  about  15  degrees.  The 
Mach  5  station  was  located  3.55  cm  from  the  throat  where  the  area  ratio  A/A*  was 
approximately  20.  The  total  length  of  the  nozzle  was  5  cm.  The  combination 


germanium  window  mounts  and  "ears"  were  machined  from  brass  and  screwed  into  the 
nozzle  flange  until  flush  with  the  sides  of  the  nozzle  blocks.  The  ears  were 
used  to  eliminate  the  effects  of  the  shock  tube  recoil  from  disturbing  the  posi¬ 
tion  of  the  laser  probe  beam  with  respect  to  the  nozzle  test  station.  The 


germanium  windows  were  one  inch  in  diameter  and  1/8  inch  thick  and  were  held 
in  position  on  the  ends  of  the  window  mounts  by  0.020  inch  thick  face  plates. 
The  windows  were  centered  *he  Mach  5  station  of  the  nozzle  and  had  a  view 
dimension  of  one  half  inch  in  diameter. 


The  laser  optical  system  was  designed  to  substantially  reduce  noise  feed¬ 
back  to  the  laser  cavity.  This  noise  feedback  was  produced  by  stray  reflections 
from  the  optical  parts  of  the  nozzle  test  station  (window  and  ears)  that  would 
vibrate  during  a  test  run.  The  noise  feedback  would  be  greatly  amplified  by  the 
laser  cavity  and  would  prevent  consistent  and  reliable  gain  measurements.  To 
substantially  reduce  the  feedback,  the  laser  cavity  was  isolated  from  the  nozzle 
optical  parts  by  inserting  a  10%  transmitting  beam  splitter  betweer  the  laser 
output  and  the  nozzle.  Thus,  only  10%  of  the  total  laser  output  of  10  watts 
was  utilized  to  make  the  gain  measurements.  The  remaining  laser  output  was 
dumped  into  a  nearly  perfect  absorber  of  10. 6y  laser  radiation  -  a  glass  of 
water.  Approximately  3%  of  the  dumped  beam  was  deflected  by  a  beam  splitter  to 
an  infrared  detector  that  was  used  to  monitor  the  laser  output  during  a  test  run. 
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The  laser  gain  measuring  beam  passed  through  the  test  station  and  into  the 
gain  measuring  detector.  A  solenoid  operated  single  sweep  chopper  was  located 
on  the  laser  side  of  the  test  station  and  was  triggered  to  chop  the  beam  after 
the  useful  test  time  was  over.  This  provided  a  measure  of  the  total  laser  out¬ 
put  that:  was  necessary  to  determine  the  gain.  The  chopper  was  necessary  to 
eliminate  uncertainty  in  the  total  laser  output  due  to  long  time  variations  of 
the  laser  power  end  effects  produced  by  loud  acoustic  and  mechanical  noise 
created  by  firing  the  shock  tube  that  would  disturb  the  laser  output  power. 

These  acoustic  and  mechanical  noises  would  typically  reach  the  laser  a  few  milli¬ 
seconds  after  useful  test  time. 

The  gain  detector  and  monitor  detector  voltage  outputs  and  the  stagnation 
pressure  were  recorded  on  oscillographs  for  each  test  run.  The  gain  detector 
output  was  recorded  on  both  a  1  msec/cm  and  2  msec/cm  time  scale  and  a  few 
jnvolt/cm  and  50  ravolt/cm  voltage  scales  respectively.  The  stagnation  pressure 
was  recorded  on  the  same  oscillograph  as  the  fast  time  scale  gain  measurement  to 
correlate  the  gain  detector  output  with  the  pressure  history.  The  monitor  detec¬ 
tor  output  was  recorded  on  the  same  oscillograph  as  the  slow  time  scale  gain 
detector  output  with  a  voltage  scale  of  5  mvolts/cm.  Monitor  detector  output 
could  then  be  compared  directly  with  the  gain  detector  output  during  a  test  run 
to  determine  the  presence  of  any  possible  feedback.  Thin  film  heat  transfer 
gauges  were  used  to  measure  the  incident  shock  Mach  number. 

A  typical  test  run  is  shown  in  Figures  3a  and  3b  for  the  5.18%  C02  gas 
mixture  at  a  stagnation  temperature  of  2300°K  and  a  pressure  of  39.8  atm.  In 
Figure  3a,  the  upper  trace  is  the  fast  time  scale  gain  trace  and  the  lower 
trace  is  the  stagnation  pressure  trace.  The  arrival  of  the  incident  shock  and 
the  sudden  jump  in  pressure  due  to  the  reflected  shock  are  clearly  visible  in 
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the  pressure  trace.  The  gain  trace  indicates  .  vary  short  rf 

«  .1  msec)  as  the  nettle  flow  gets  started  and  then  the  production  of  positive 

gain  as  the  flow  Is  established  at  the  nettle  Mach  5  station.  The  pressure  drops 

%  msec  after  the  reflected  shock  wave  pressure  rise  due  to  the  arrival  of  the 

Shock  tube  driver  expansion  fan.  The  gain  correspondingly  goes  down  to  tero  or 
absorption. 

In  Figure  3b  the  upper  trace  Is  the  slow  tin*  scale  gain  trace  illustrating 
the  gain  signal  described  above  and  the  chopping  of  the  CO.,  gas  laser  probe  signal 
by  the  single  sweep  chopper  to  determine  the  relative  power  of  the  laser  probe 
signal  With  respect  to  the  gain  signal.  The  lower  trace  is  the  monitor  detec¬ 
tor  and  it  indicates  that  virtually  no  feedback  or  noise  affected  the  probe 

laser  during  the  test  run.  The  measured  gain  per  meter  for  this  particular  test 
run  was  0.55/meter. 


Results  and  Discussion 

The  measured  gain  vs.  stagnation  temperature  is  shown  in  F1gure  4.  The 
temperature  for  each  data  point  was  calculated  on  the  basis  of  the  measured  in¬ 
cident  shock  Mach  number  and  equilibrium  flow  theory  for  a  reflected  shockwave 
without  chemical  reactions.  The  theoretical  curves  are  based  on  the  kinetic 
model  and  rate  data  as  described  in  References  5  and  6.  This  kinetic  model  and 
rate  data  have  been  shown  to  be  valid  by  experimental  verification  and  comparison 
other  works  (References  4  and  7)  at  temperatures  below  2000°K. 

The  experimentally  measured  gain  data  agree  with  theory  at  2000°K.  At  tem¬ 
peratures  above  2000%  the  measured  gain  deviates  considerably  from  the  theory 
for  the  2.63%  aad  5.18%  CO,  gas  mixtures  and  only  slightly  for  the  10.73%  CO, 
mixture.  A  maximum  in  the  measured  gain  occurs  at  about  2500°K  for  both  the' 
2.63%  and  5.18%  CO,  gas  mixtures,  the  gain  at  these  maxima  being  about  50%  and 
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and  30%  higher,  respectively,  than  the  theoretically  predicted  gain.  Above 
2500°K  the  measured  gain  gradually  decreases  but  still  remains  considerably 
higher  than  theory.  One  test  run  (not  shown  on  Figure  4)  was  performed  with  the 
2.61%  C02  gas  mixture  at  a  stagnation  temperature  of  5000°K  which  indicated  zero 
gain.  Possible  mechanisms  that  could  account  for  this  difference  between  theory 
and  experiment  are  discussed  below. 


Figure  5  is  a  plot  of  gain  divided  by  CO^  concentration  vs.  stagnation 
temperature.  This  plot  shows  that  the  principle  of  obtaining  higher  gain  per 
C02. particle  by  using  starved  C02  gas  mixtures  is  applicable  for  high  tempera¬ 
tures  as  well  as  high  pressures.  At  2500°K  the  gain  per  particle  for  the  2.63% 
C02  mixture  is  over  4  times  larger  than  the  gain  for  the  10.73%  mixture  and  1.4 
times  larger  than  the  5.18%  mixture.  This  increase  in  gain  per  C02  particle  is 
brought  about  by  the  increased  efficiency  of  the  freezing  process.  The  effi¬ 
ciency  of  the  freezing  process  is  characterized  by  the  equation  for  gain. 


g 


NC02 

N 

TOTAL 


(Xu  “  V 


where  g  gain,  number  density  of  C02>  =  total  number  density  of 

gas  mixture,  and  ^  =  fraction  of  C02  molecules  in  the  upper  and  lower 

laser  levels,  respectively.  We  can  write  the  mole  fraction  of  C02  in  the  gas 
mixture,  as  \ PCq^  =  NC02^NT0TAL  an^  hence  obtain 


Thus  we  see  that  g/V'CQ2  is  a  more  direct  measure  of  the  population  inversion 
than  g  is.  The  larger  g/i| >co2  is,  the  larger  the  population  inversion  is,  hence 
the  more  efficient  the  freezing  process  is.  Energy  storage  of  the  gas  mixture 
is  also  increased  since  a  larger  fraction  of  the  C09  and  N2  is  in  the  upper 
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excited  state. 


It  is  apparent  from  the  plot  of  the  data  in  Figure  4  that  the  mechanism 
that  io  producing  the  discrepancy  between  theory  and  experiment  is  very  effec- 
tive  on  the  2,63%  311(1  5.18%  C02  mixture  but  is  rather  ineffective  on  the  10.73% 
C02  gas  mixture.  Thus,  the  mechanism  must  display  not  only  a  temperature  sensi¬ 
tivity,  but  also  a  very  strong  sensitivity  to  the  concentration  of  C02.  This 
sensitivity  of  the  results  to  the  concentration  of  CO,  will  prove  to  be  the 

major  stumbling  block  to  any  "simple"  explanation  of  the  discrepancy  between 
theory  and  experiment . 

Possible  uncertainties  in  the  rate  data  for  the  reaction 
C02  +  N*  $  C02(v3)  +  N2 

were  looked  into  as  a  possible  source  of  error  in  the  theoretical  predictions. 

_ate  data  for  this  reaction  is  not  accurately  knovm  for  temperatures  near 
o  (5) 

2000°K.  The  rates  were  varied  by  as  much  *s  factors  of  10  at  high  tempera¬ 
tures,  but  no  consistent  and  reasonable  results  were  obtained  that  matched  the 
experimental  data  points. 

The  possible  effects  of  the  upper  excited  vibrational  energy  states  of  N2 
adding  to  the  first  vibrational  energy  state  can  be  estimated  on  the  basis  of  a 
simple  calculation.  Assuming  equilibrium  at  the  throat,  we  can  calculate  the 
number  densities  of  the  N2  molecules  in  the  various  vibrational  energy  states 
at  the  throat..  If  we  assume  that  all  the  nitrogen  in  the  v  =  2  excited  state 
cascades  down  to  the  v  -  1  state  and  remains  frozen  as  the  flow  expands  through 
the  nozzle,  we  can  recalculate  the  vibrational  temperature  of  the  N2  molecules 
by  simply  adding  the  number  density  of  these  extra  excited  molecules  to  the 
number  of  v  =  1  exc. ted  molecules  as  calculated  from  the  regular  theory.  The 
following  tat  e  shows  the  results  for  a  stagnation  temperature  of  2500°K: 
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Gas  Mix 


10.73  CO, 


5.18  CO, 


2.65  CO, 


mlar  Theoi 


Tvib  Gain 


1100  .36 

1500  .45 

1750  .31 


Simple  Addition  Experimental 


Process 


Gain 


We  see  that  the  correction  is  reasonable  for  the  5.18?  CO.,  mix,  a  little  low  for 
the  2.632  C02  nix,  but  completely  wrong  for  the  10.73*  mix.  The  10.73*  „,ix  may 
deactivate  the  upper  excited  vibrational  states  faster  which  may  make  this  simple 
calculation  Loo  high  a  correction.  However,  the  calculation  for  the  2.632  CO., 
gas  mixture  would  seem  to  indicate  that  the  simple  cascade  process  would  have^o 
be  1002  efficient,  a  very  unlikely  result.  Also  this  simple  mechanism  does  not 
provide  for  a  maximum  gain  at  about  2500»K.  It  would  seem  unlikely  that  this 
process  can  account  for  the  observed  discrepancy  between  theory  and  experiment. 

Dissociation  effects  can  manifest  themselves  in  three  ways:  (1)  Chemistry 

effects  behind  the  reflected  shock  wave  will  reduce  the  stagnation  temperature 

for  a  given  initial  pressure  and  incident  mach  number,  (2)  dissociation  of  CO,, 

may  self  starve  the  gas  mixture  during  the  early  portion  of  the  freezing  process, 

(3)  the  energy  of  recombination  may  be  dumped  into  the  excited  vibrational  energy 
state  of 

To  determine  the  effects  of  CO.,  dissociation  on  the  stagnation  temperature, 
chemical  equilibrium  calculations  for  a  reflected  shock  wave  were  performed  as¬ 
suming  the  following  seven  chemical  species  were  present,  CO.,,  N  He,  CO,  0  ,  0, 
and  HO.  The  chemical  equilibrium  temperature  behind  the  reflected  shock  wave  was 
calculated  lor  a  given  pl>  and  and  compared  to  the  no-chemistry  calculations. 
The  following  table  shows  the  results. 
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/ 
t . 


Gas  Mix 


10.73%  CO, 


T0  (no  dissociation) 


2000 

2500 

3000 


To  (dissociation)  degree  of  dissociation 


1980 

2430 

2750 


.008 

.063 

.180 


5.18%  CO, 


2000 

2500 

3000 

3500 


1990 

2430 

2830 

3160 


.016 

.090 

.280 
.  ^80 


2.63%  CO, 


2000 

2500 

3000 

3500 


2000 

2490 

2890 

3270 


.018 

.140 

.390 

.660 


These  results  Indicate  no  significant  difference  in  the  stagnation  temperature 

•  20°°°K  "here  tha  deerea  °£  dissociation  of  N.,  is  very  small.  Above  1000%  the 
T0  (dissociation)  is  lower  than  the  Tq  („„  dissociation),  ..he  difference  amounting 
to  about  340°K  at  35001  far  the  5.18%  CO,  mix.  It  should  also  be  noted  that  the 
error  in  determining  To  is  on  the  some  order  as  the  temperature  correction.  Cor¬ 
recting  the  experimental  data  on  Figure  3  for  the  lower  stagnation  temperature 
would  nicely  correct  the  experimental  results  for  the  10.73%  CO.,  mix.  However, 

U'  IS  oWl0u!i  that  lmrarl"8  the  temperature  for  the  2.63%  and  5.18%  CO.,  gas  mix¬ 
tures  would  not  correot  the  data  (except  at  the  highest  temperatures  -  3500)  as 
the  measured  gain  between  2000”K  and  30C0"K  lie  entirely  above  the  maximum  theo¬ 
retically  predicted  gain.  Some  other  mechanism  is  affecting  the  gain  for  the  two 
lower  C02  gas  mixtures.  The  effects  of  self  starving  <s  not  expected  to  be  of 
any  importance  due  to  the  low  degree  of  dissociation  of  the  CO.,,  particularly  at 
2500°K  where  there  is  a  maximum  in  the  measured  gain.  If  chemical  equilibrium 
is  assumed  at  the  throat  of  the  noxxle,  the  degree  of  dissociation  drops  by  about 
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I 


a  tact->r  of  4  and  hence  is  even  more  Insignificant  as  far  as  the  self  starving 
effect  is  concerned.  '  . 

!  I  1 

The  amount  of  energy  released  in  the  recombination  of  CO  and  0  to  form 

'  :  I  t  'I  ’ 

C02  would  have  a  significant  effect  on  the  gain  if  it  were  dumped  into  the  first 

1  '  1 

vibrational  energy  level  of  N  ,  even  though  the  degree  of  dissociation  is  very 

i  ,  ' 

small.  In  order  for  the  recombination  process  to  add  energy  to  the  vibrational 

i  i  'i 

energy  levels  of  N2  in  such  a  manner  that  the  energy  is  frozen  in  these  vibra- 

'■  i  i  - 

tional  energy  states,  the  relaxation  time  fof  recombination  must  be  on  the  same 

5  •  t  » 

order  as-  or  shorter  than  the  vibrational  relaxation  time  for  N  .  Unfortunately, 
the  following  calculation  will  indicate  chat  the  recombination  process  is  essen¬ 
tially  frozen  throughout  the  supersonic -.portion  of  the  nozzle  and  hence  cannot  ‘ 
contribute  any  energy  to  the  vibrational  energy  states. 

Davis  has  reported  that  the  forward  rate  constant  fop:  the  reaction 

kf  !  ,  \  i  I 

c°2  +  N2  CO  +  0  +  N2  1  1 

:kb  '  '  '• 


G 


-  4.07  x  10  exp(0  J. 23eV/kT)  cm3particle  *sec  11 


is  given  by 
ki 

The  backward  rate  can  be  obtained  by,  using  the  equilibrium  constant  for  the  above 

I 

*  «  !  I 

reaction.  For  the  5.18%  C02  gas  mix  at  a  stagnation  '‘empnrature  of  2500°K  at, 40 
atm,  the  equilibrium  temperature  at  the  'throat  is  2100 °K  at  21  atm.  Chemical 
equilibrium  calculations  gives  u  ,  the,  following  number  densities  at  the  thtoat: 


N._  =  3.5  x  10la 
N2 

nco2  =  3-7  x  1018 


NCq  =  1.04  x  10 
Nq  =  5.0  x  10 llf 


1  7 


nno  r  2-a  51 1017 


,NHe  -  5,5  x  10 


19 


N  -  3.7  x.  1016 
U2 
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Tilt  relaxation  time  for  the  recombination  process  is  given  by 

TrecK  V^=2'3mSeC 

An  upper  limit  for  the  N2  vibrational  relaxation  process  can  be  obtained  by  con 
sidering  a  pure  N2  gas  mix.  ?rom  Millikan  and  White ^  we  obtain  for  2100°K, 
pTyib  »  500  atm  ysec, 


hence 

l 

.  i  __  500 

!  *Tvib  *  IF  "  2^sec 

Henc-  Trec  >>  TvjQ,  an^  is  unlikely  that  recombination  in  the  supersonic  portion 
of  the  nozzle  will  have  any  affect  on  the  vibrational  temperature  of  nitrogen. 

Similar  results  were  obtained  for  the  other  C02  gas  mixtures. 

! 

i  If  the  ener«y  of  recombination  is  dumped  into  the  vibrational  energy  level 
sub sonic ally  (the  degree  of  dissociation  varies  from  about  10%  to  3%  between  the 
plenum  and  the  throat)  the  vibrational  temperature  of  N2  will  be  higher  at  the 
throat  than  the  transit,  .ional  temperature. 

Calculation:;  were  performed  using  the  regular  theory  to  determine  the  effect 
on  the  gain  of  setting  the  N2  vibrational  temperature  at  the  throat  500°K  higher 
than  the  equilibrium  temperature  of  2100°K.  These  calculations  indicated  that 
i  the  'excess  vibrational  energy  was  quickly  dumped  back  into  translational  energy 
and  lower  C0?  vibrational  energy  levels  by  the  relaxation  collisions  between 
C02  and  N2*  The  net  result  was  virtually  no  change  in  the  predicted  gain.  Thus 
it  is  concluded  that  on  the  basis  of  the  above  precursory  analysis  dissociation 
has  a!  negligible  effect  on  the  vibrational  relaxation  process  and  gain. 

The  contributions  of  other  near  resonance  transitions  may  contribute  sig¬ 
nificant!  v  to  the  predicted  gain.  The  predicted  gain  was  calculated  upon  the 
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basis  of  the  C02  vibration  -  vibration  transition  001  ->  100.  Other  near  resonant 
transitions  can  be  represented  in  general  by 
m,  n,  p+1  -*■  m+1,  n,  p 

where  m  =  0,  1,  2,...,  p  =  0,  1,  2,...,  and  n*=  0,  1,  2,  3 .  If  we  assume 

that  these  transitions  are  exactly  resonant  with  the  001  100  transition  we  can 

calculate  the  maximum  possible  effect  on  the  predicted  gain. 

The  predicted  gain  for  a  given  transition  is  proportional  to  the  following: 


cc _ l _  Nm,n,p+1  Nmfl,n.p  I 

a,n,p  Vp+iL  n  n  J 


where  N  *=  total  number  of  C02  molecules, 

Nra,n,p+1  *  number  of  CO,,  molecules  in  the  m, n, p+1  state  or  upper 
laser  level. 

Nm+l,n,p  "  n,1T7lber  of  molecules  in  the  lower  laser  level, 
fcm,p+l  “  lifetime  of  the  upper  la-er  level. 

Based  upon  harmo.de  oscillator  theory,  we  assume 

t  cc  —  1 

m,p+l  (m+1)  (p+1) 

The  total  gain  due  to  all  possible  transitions  is  given  by 

_  99  99  oo 

G  »  XL  JL  z  0 

m=U  n=U  p=0  °m,n,p 

Basov  et  ai:9  gives  for  the  quantity  N  /N  the  following  expression 


Nm,n,p  r  -m0!/1,  “0i/Ti  -nG^/T,  -0,/T„  ~p0  /T  -0  /T 

~N  "  Ce  (1-e  )][(n+l)e  2(l-e  2)2]fe  3  3(l~e  3  3)] 


where  0.^  02,  03  are  the  characteristic  temperatures  for  the  o  ,  v»,  V  vibrationa: 

**  *>  J 

modes  of  C02  <1920°K,  960°K,  and  3380  K  respectively),  and  T. ,  T„,  T  are  the 

~  J 

effective  vibrational  temperatures  for  the  modes  V^,  v2>  v3  respectively. 
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Putting  Basov;s  expression  into  the  expression  for  the  total  gain  and  summing  we 
obtain 


■e,/x  -62/t2 

)  (1-e 


-ejT 

)  (l-e  3)]'2 


where  gQ  -  gain  for  the  001  100  transition.  This  correction  factor  to  the  pre¬ 

dicted  gain  gQ  is  illustrated  in  figure  6  for  the  2.63%  C02  mixture.  The  tempera¬ 
tures  T1  and  T3  are  the  temperatures  that  the  upper  and  lower  laser  levels  (001 
and  100  respectively)  froze  at  as  determined  from  the  theory  for  predicting  g  The 
correction  to  the  experimental  data  points  from  the  lowering  of  the  stagnation 
temperature  due  to  dissociation  effects  are  also  illustrated.  On  the  basis  of  a 

4  temperatur-  model  for  the  kinetics  of  a  C0„-N„-He  gas  mixture  we  assume  T  **  T  . 

£  c  12 

The  correction  factor  to  the  gain  g^  U  *bout  a  factor  of  two  larger  than 
the  experimentally  measured  gain.  However,  the  correction  factor  does  indicate 


a  max  ..mum  at  about  250C°IC  and  the  shape  of  the  curve  is  in  fair  agreement  with 
the  experimental  data.  (This  can  be  easily  seen  by  lowering  the  curve  for  G  such 


that  G  =  gQ  at  200G°K) . 


Similar  results  were  obtained  for  the  5.18%  CO,,  gas  mixture.  Results  for 
the  10.73%  gas  mixture  indicated  no  maximum  in  the  gain  and  the  correction  factor 
gave  results  that  were  only  about  1.3  times  larger  than  the  experimental  data. 

This  latter  results  was  primarily  due  to  the  smaller  differences  between  ^  and  T3 
due  to  the  faster  relaxations  effects  for  the  10.73%  C02  mix.  All  of  these  results 
are  in  good  qualitative  agreement  with  the  experimental  data. 


That  the  correction  factor  gives  a  gain  that  is  too  large  is  not  surprising 
siAce  exact  resonance  was  assumed.  The  magnitude  of  the  correction  factor  would 
be  substantially  reduced  due  to  the  fact  that  the  various  trar  Itions  are  not 
exactly  resonance.  However,  the  general  shape  of  the  curve  is  expected  to  stay 
the  same,  i.e.  showing  a  maximum  at  about  2500°K.  Calculations  performed  by 
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tlunj ee  and  Christiansen  (see  appendix  A)  on  the  effects  -f  pressure  broadening  of 
these  near  resonance  lines  on  the  absorption  coefficient  at  10.6,1  indicates  that 
the  effect  of  these  near  resonance  lines  is  typically  1/3  as  effective  a  exact 
resonance  could  be.  The  results  of  reducing  the  difference  between  G  and  gQ  by 
a  "actor  of  1/3  (i.e.  gain  -  g()  +  1/3(G  -  g^)  is  aIs0  lllustratcd  la  flgJe'  6. 

Aj. though  good  agreement  is  obtained  for  T,  >2500%  the  predicted  gain  is  still  too 
high  for  Tq  <  2500°K.  A  maximum  at  about  2500°K  is  still  indicated. 

Qualitatively  this  near  resonance  effect  appears  to  offer  a  very  satisfac¬ 
tory  explanation  for  the  difference  between  the  measured  and  predicted  gain. 

More  detailed  calculations  are  in  the  process  of  being  formulated  and  a  experi¬ 
mental  program  involving  changes  in  pressure  as  well  as  temperature  is  being 
net  up  to  more  firmly  establish  (or  disprove)  that  this  near  resonance  effect  is 
'-.deed  t;ie  mechanism  that  is  responsible  for  the  discrepancy  between  predicted 

and  measure  gain.  The  results  of  these  investigations  will  be  reported  in  the 
near  future. 


III.  ENTROPY  PRODUCTION  DUE  TO  RADIATION  EXTRACTION 

Preliminary  analysis  '  of  the  thermodynamics  of  lasing  a  flowing  media 
has  shown  that  the  normal  dissioative  effects  and  entropy  production  can  be  re¬ 
duced  ideally  to  zero.  This  analysis  wa  based  on  a  two  temperature  model  at 
low  concentration  of  C02  in  which  the  up* or  la  er  level  of  C02  and  the  vibrational 
energy  of  could  be  characterized  by  a  single  temperature,  T  .  The  lower  laser 
level,  the  lower  vibrational  states,  and  translational-rotational  heat  bath  were 

characterized  by  another  temperature  T.  For  a  one  dimensional  flowing  gas  with 
no  collisional  deactivation,  it  is  easy  to  show  that:  a  small  variation  in  entropy 
of  the  fluid  during  lasing  is  given 
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6S  «.  -  l]  ~  »  [l  -  (i~n)~]6sv 


where  n  is  the  quantum  efficiency  for  C02(0.4)  and  6s  and  6sv  are  the  variation 
in  entropy  of  the  system  and  vibrational  states  respectively,  6q_  is  the  radiated 
energy.  The  requirement  of  gain  in  the  .laser  system  also  stipulates  that 


T  /T  >  ~~  -  1. 
v'  -  1-n 


68 


Thus  for  all  T^/T  >  it  is  seen  that  6s  >  0,  since  6s^  <  0  when  lasing  occurs 

(for  Tv/T  =  or  zero  gain,  6s  -  0) .  For  Ty/T  slightly  greater  than  6s 
is  only  slightly  greater  than  0  and  it  can  be  shown  that 


6s  = 

ue 
_v 

v 


6e  2 
v 


6V 


cT  2  cn2T  2 
v  v 


where  c  -■  dT  and  6ey  is  the  variation  in  vibrational  energy. of  N„.  In  this 

V  «- 

case  entrop:  production  is  a  higher  order  phenomena  than  laser  energy  out.  Thus 
laser  cavity  systems  with  finite  power  extraction  and  no  entropy  production  can 
be  designed  in  principle.  That  is 

?6c2 


AS  ..  Pi  1  Pli 

^  P%±-  ^ 


0 


^Li  rMLi 

where  AS  and  AQ^  are  the  net  entropy  change  in  the  fluid  and  net  finite  laser 


energy  output.  To  show  this,  let  6q  -  6q  =  ^ 


m 


as  a  constant  differential 


energy  output.  Then  the  above  equation  becomes 


lim  AS 

m  ^  oo  — - 


Vm  AQ 

~7  0 


AQj  cT26q^m  cT2m 

In  the  limit  of  Infinite  number  of  steps,  entropy  production  goes  to  zero  with 
finite  AQ  . 

Lt 
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Io  achieve  this  ideal  result  in  prentice  requires  that  lasing  takes  place 

in  a  variable  area  channel  such  that  Ty/T  remains  at  all  times  close  to  the  zero 

gair /meter  value.  For  the  case  of  T/T  equal  to  a  constant  we  have  shown  that 

the  entropy  of  the  fluid  is  related  to  the  entropy  change. of  the  vibrational 
states  by 

AS  =  [l  -  (l-n)T  /T]|  AS 
v  I  v 

•  1  const  T  /T 

v 

Again  we  see  that  if  Ty/T  =  =  const.,  AS  =  0. 

There  are  a  number  of  questions  to  be  asked  concerning  the  validity  of  the 
simplified  result  for  a  CO,,  mixture  before  detailed  examination  of  a  cavity 
design  can  begin.  Firstly,  is  mode  analysis  a  valid  approximation  for  C02?  The 
general  equations  governing  the  transfer  of  energy  between  the  vibrational  mode 
of  N2,  the  three  modes  of  C02 ,  and  translational-rotational  bath  are  similar  to 
those  used  in  reference  9  with  the  addition  of  radiation  terms.  However,  at  high 
radiation  intensities  the  energy  transfer  between  modes  is  much  more  rapid  than 
that  due  solely  to  collisional  processes.  Furthermore,  radiation  transfer  usually 
occurs  between  two  particular  vibrational  levels  and  not  between  modes.  The 
validity  of  a  mode  analysis  rests  on  the  basic  assumption  that  intramode  V  -  V 
transitions  are  sufficiently  rapid  to  ensure  a  Boltzmann  distribution  of  levels 
at  a  characteristic  mode  temperature  within  each  mode.  This  may  not  be  valid 
at  high  radiation  intensities  and  very  low  C02  concentrations.  However  while 
these  probabilities  are  not  precisely  known  it  is  believed  that  V  -  V  energy 
transfer  probabilities  may  be  very  large  0(%)(ll>  and  therefore  allow  the  use  of 
mode  analysis  down  to  a  level  of  CC>2  concentration  of  about  0.1%,  Furthermore, 
some  analysis  by  Yesil  (Appendix  B)  shows  that  even  if  V  -  V  rates  are  small  or 
even  ignored,  it  is  difficult  to  get  a  very  non-Boltzmann  distribution  of  vibra- 
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..lonal  states  for  time  scales  typical  of  GDI's.  Therefore  we  feel  mode  analysis 
is  a  valid  model, 

*  «•  .  •  . 

Another  important  question  to  ask  is  how  does  the  two  temperature  model 
compare  with  a  more  realistic  4-temperature  model  of  Munjee^  and  others  (see 

also  Appendix  C)  (T  ,  T  T  =  T  -n  r  *.1. 

vN„  v3C0  ’  v  1v  *  0ne  of  the  concerns  is  the  ef- 

—  2  12 

feet  01  radiative  saturation  on  entropy  production.  Unless  intensities  near 
saturation  (defined  as  intensity  required  to  reduce  gain  by  two  in  presence  of 
radiation)  are  used,  the  laser  medium  will  be  ineffective  and  collisional  decay 
will  result  in  a  loss  of  vibrational  energy  and  consequent  increase  of  entropy. 

if  high  intensity  radiation  is  introduced  to  the  system,  then  the  mode 
temperatures  are  affected,  also  causing  entropy  production.  Part  of  this  pro¬ 
blem  also  can  be  estimated  by  the  following.  The  entropy  of  the  fluid  for 
a  4-temperature  model  is 


6s 


6e 


de  dp 

TR  vCO„  vCO„  de 

~  +  Pdv  +  +  - 1  + 

T2  T3  TvN 


where  devC0^/T2  ie  the  contribution  due  to  modes  1  and  2  of  C02, 


de 


vCO, 


the 


de 


vN 


entropy  contribution  due  to  code  3  of  CO  and  Ja  the  contribution  of  N„.  If  „e 
specialize  to  low  concentration  of  C02,  then  the  C02  contribution  to  the  entropy 
of  the  system  can  be  neglected,  hence 

de  +  pdv  de 

6s  *  — ±£ _ _ _ $1 

T  T  . . 

vN 

From  reference  10,  this  is  still  the  same  formula  for  the  two  temperature  system 
and  thus 


X  d 

4»  "  hi-n>-^  -  i]  ^  -  Ei  -  <i-n>^  ]  6sv 


N 
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Ho»  in  this  case  T  <  ^  and  >  Tj  (see  e.g.  10)  as  the  tore  gain  condition  is 
VT2  ’h'  1'68'  Therefore  6s  >  0  under  all  cases  and  a  net  entropy  is  pro¬ 
duced.  We  can  readily  see  with  the  help  of  reference  12  that  as  a  laser  mixture 
is  radiatively  saturated  I  «  r2  and  »  l3  for  large  gain  cases  that  a  large 

entropy  is  produced.  However  for  a  low  gain  case  I  and  will  be  closer  to  T, 
and  ’t'3  respectively  and  such  large  deviations  will  not  be  possible.  Now  the 
rate  of  power  extraction  is  given  by 
dl/dx  =  gl  (power/voluiue) 

In  the  limit  of  very  high  intensities,  the  rate  of  power  extracted  It  a  maximum 
(for  homogeneously  broadened  gain)  and  given  by 

a  T 

dx  I  ■*  oo  g0g 

Since  Ig  does  not  depend  on  the  vibrational  energy  of  the  systen/12)  dl/dx| 

I  ^  CO 

is  a  smaller  rate  for  the  low  gain  case  (because  being  deliberately  kept 

near  0)  than  the  high  gain  case.  Thus  in  the  presence  of  collisional  deactiva¬ 
tion  of  the  laser  states  we  have  a  compromise  situation  developing.  By  reducing 
the  unsaturated  gain  toward  aero  we  reduce  the  entropy  production  due  to  radia¬ 
tion  extraction,  but  the  rate  of  extraction  is  also  reduced  so  that  collisional 
decay  becomes  important.  Thus  the  optimum  in  saturation  is  one  which  reduces 
gain  only  to  a  point  where  the  net  entropy  generation  is  the  smallest.  Estimates 
show  that  the  minimum  value  of  g^  is  proportional  to 


i  (St.  +  t. 
'min  -  — - - i 


4  fn(p) 


where  the  coefficients  axe  defined  in  reference  12,  The  bajic  relaxation  times(6) 
are  sensitive  functions  of  temperature  and  larger  values  of  gQ  will  be  required 
at  higher  translational  temperatures. 
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Some  preliminary  calculations  using  the  computer  simulation  for  various 


A' 


nozzle  geometries  are  shown  in  the  table  below.  In  all  the  cases  below,  the 


medium  was  fully  saturated,  i.e.  gain  was  kept  at  zero  where  lasing  occured. 


EXAMPLE  OF  ENTROPY  PRODUCTION  DURING  LASING 


Basic  Conditions 
sntratio 
100  atm 


Concentrations:  1%  C02>  60%  N2,  39%  He 


T  =  1500°K 
o 


I  for  g  =  0 
o 


Throat  height,  h*  =  0.1  cm 


CASE  1 


Nozzle  Description 
20°  half  angle  wedge  to  A/ A*  =  10, 
Lase  at  constant  area  channel 


CASE  2 


20°  half  angle  wedge  to  g^  =  0 


10°  to  A/A*  =4.1 
5°  to  A/A*  «  6.2 
2°  to  A/A*  -  23 


CASE  3  20°  half  angle  wedge  to  g  =0 


6°  to  A/A*  =3.6 
3°  to  A/A*  =  5.1 
1°  to  A/A*  =  13 


Ae 


_ v 

12.9  J/gm 


12.7  J/gm 


AS/R 

.079 


.068 


12.9  J/gm  .052 


$  1 


Case  1  represents  a  standard  case  whereby  lasing  takes  place  in  a  constant  area 
channel  where  the  Mach  number  is  about  4.  During  the  lasing  phase  there  is  about 
an  8%  increase  in  entropy  for  an  energy  extraction  of  5.9  KJ /#  (12.9  J/gm).  How¬ 
ever  by  contouring  the  nozzle  differently  during  the  pox^er  extraction  phase,  as 
in  Cases  2  and  3,  an  improvement  is  possible  for  the  sane  energy  out /gram.  The 
best  case  is  number  3  where  there  is  only  about  5%  of  an  entropy  increase.  This 

is  not  as  good  as  we  would  like  and  it  certainly  does  not  represent  an  optimum 

| 

situation.  The  results  do  reflect  the  effect  of  collisional  deactivation  of  the 
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Up.  *  laser  level  and  the  4-temperature  model.  He  are  looking  Into  numerical 
procedures  whereby  the  optimum  nozzle  geometry  can  be  calculated  by  the  computer, 

V* 

IV.  NEW  NOZZLES  AND  NOZZLE-DIFFUSER  COMBINATIONS 

The  Supersonic  Nozzle- 

A  may  of  accompl..  ,hlng  the  rapid  ‘ncrease  In  flow  cross-sectional  area 
required  for  an  all  supersonic  gasdynamic  laser  Is  shown  schematically  on 
«Sure  7.  Here  an  Incoming  supersonic  stream  Is  turned  and  expanded  through  an 
array  of  simple  Prandtl  -  Meyer  expansion  fans.  One  side  of  each  nozzle  element 

consists  of  a  sharp  comer,  while  the  other  is  contoured  to  cancel  the  incident 

expansion  wave.  Due  to  the  fact  that  Qii  n,  ,  ,  . 

1  of  the  high  velocity  in  coruing  flow  is 

directly  expanded,  a  high  mass  flux  and  power  are  expected  for  a  given  pressure 

capability.  This  feature  makes  the  supersonic  turning  nozzle  particularly 
attractive  for  shock  tube  applications. ^13^ 

The  kinetic  performance  of  these  nozzles 

uesc  nozzles  has  been  compared  with  conventional 

convergent-divergent  nozzles  a  - 

8  nozzles.  A  two  temperature  model  was  used  for  the  N  -CO 

system,  and  the  calculations  were  carried  along  separate  streamlines  in  orders  ’ 
find  average  exit  properties.  It  was  found  that  an  individual  supersonic  nozzle 
will  produce  less  power  per  unit  of  mass  flux,  (X),  for  the  same  mixture,  expan¬ 
sion  area  ratio,  total  temperature  (T,) ,  and  product  of  total  pressure  (P  ) 
times  throat  height  (h, .  The  difference  le  primarily  due  to  the  lower  Jering 
vibrational  temperature  for  the  high  speed  case.  X  Is  given  by 
X  -  0.4[ev(Tv)  -  e  (1.68T)] 

2 

he  vibiational  'energy,  ey,  Is  evaluated  for  a  particular  N.,  concentration, 

Tv  rs  the  vibrational  temperature,  and  T  Is  the  translational-rotational  tempera- 

*^ure  of  the  expanded  flow.  T  anri  t  „„„  i 

v  and  T  are  known  for  perfect  freezing  (p  h  =  0) , 

eValUated  f°r  th£  nOZZle  configuration  being  consldere 
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*  .JM  ’ 


4-.;, 


as  a  first  approximation  the  ratio  of  this  parameter  at  ,  h  -  0  can  be  used  to 


p:  edict  relative  performance  at  finite  >#h.  Mete  precise  calculation  requires 


O  *  - v«*vuj.wv*wu  i.CL[  U_L  L  Oj 

the  use  of  the  two  temperature  model  calculations , (l4)  or  even  better,  the  four 
temperature  program  of  Appendix  C.  The  X  reduction  for  an  individual  supersonic 
nozzle  could  be  easily  overcome  by  a  decrease  In  relative  poh  (the' supersonic 
nozzles  operate  In  a  lower  pressure  environment  and  are  stronger),  an  Increase 


In  the  expansion  area  ratio,  or  by  a  direct  increase  in  To.  As  will  be  dis¬ 


cussed,  the  supersonic  nozzles  appear  to  be  more  suitable  for  staging  purposes 
because  of  the  relatively  low  p<>  losses.  Further,  the  power  from  a  single 
stage  Is  considerably  greater  for  the  turning  nozzle  because  the  number  of  noz¬ 
zle  elements  that  can  be  fitted  onto  a  given  duct  is  increased  directly  as  the 
nozzle  area  ratio. 

A  grid  of  supersonic  turning  nozzles  has  been  constructed,  consisting  of 
twenty  h  -  1  mm  elements  which  together  expand  a  2  x  5  cm  flow  to  JO  x  5-cm. 


The  array  Is  fed  through  a  sharp-edged  rectangular  pipe  by  a  7.5  cm  shock  tube. 

(15) 


Using  an  external  C02  laser,  gain  has  now  been  measured  in  the  expanded  flow 
It  has  been  found  necessary  to  shorten  the  inlet  pipe  and  use  the  two  extreme 
nozzles  as  boundary  layer  bleeds.  This  stopped  choking  of  the  end  nozzles 
and  consequent  propagation  of  shock  waves  across  the  cavity  flow,  a  sample  gain 
trace  Is  shown  on  figure  7  for  a  run  with  the  boundary  layer  bleeds  Installed. 


A  mixture  of  40%  Nj,  10%  C02,  and  50%  He  by  volume  was  used,  with  the  , shock  tube 


operated  so  as  to  produce  To  -  2100%  P<)h  «  1.0  atm  cm,  and  an  inlet  Mach  numbe, 


Of  *  1.5.  Absorption  is  first  observed  as  the  incident  shock  begins  to  cut  acrot 
the  beam  axis.  From  earlier  aerodynamic  studies (16>,  „  rearward  facing  starting 
Shock  is  expected  to  follow  and  Initiate  the  lasing  flow.  This  wave  begins  to 
cross  the  beam  axis  after  .  200psec  and  the  absorption  decreases  as  the  noneqni- 
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librium  gas  appears.  The  lower  trace  on  the  figure  was  obtained  with  a  pressure 
gage  mounted  in  the  cavity  wall  and  more  clearly  shows  the  passage  of  the  two 

1r  i?}  % 

shock  waves.  The  second  or  starting  shock  is  slow  to  form  and  move  out  of  the 
cavity,  and  little  if  any  plateau  is  reached  in  the  gain  trace  before  the  flow 
begins  to  be  cut  off  by  the  arrival  of  the  shock  tube  driver  gas.  The  gain  climbs 
near  to  the  calculated  level,  but  the  experiment  is  not  definitive.  As  calcula¬ 
tions  show  that  lengthening  the  shock  tube  will  not  substantially  increase  the 
test  time,  a  conversion  to  combusted  Ludwieg  tube  operation  is  being  carried  out. 


It  was  mentioned  in  the  Introduction  that  nozzle  arrays  can  be  staged  so  as 
to  increase  X  if  sufficient  pressure  ratio  is  available.  As  a  first  approxima¬ 
tion  to  a  staged  system,  assume  that  pQh  is  kept  constant  for  each  stage  by 

■» 

increasing  the  nozzle  dimension  as  pQ  falls.  Assume  also  that  the  laser  energy 
extracted  is  a  small  fraction  of  the  total  enthalpy  so  that  Tq  is  essentially  ex¬ 
changed.  If  A  is  the  pQ  entering  a  stage  divided  by  that  recovered  at  the 
diffuser  exit,  then 

An  =  p  /p 

roi  roe 

where  n  is  the  number  of  stages  that  can  be  accommodated  between  the  total  pres¬ 
sure  entering  the  first  stage  (Poj),  and  that  exiting  the  last  stage  (p  ) .  The 

specific  power  of  a  staged  system,  X  is  then  given  by  X  times  n  or 

s 


Xs  =  (X/lnX)  In  Poi/poR 


where  X  is  the  specific  output  per  stage.  The  performance  of  a  staged  system  thus 
depends  on  PQi/poe  and  X/lnA,  the  latter  being  a  staging  figure  of  merit. 

It  was  just  pointed  out  that  X  is  reduced  for  the  all  supersonic  nozzle. 

The  nonequilibrium  pQ  loss  would  also  be  reduced,  but  more  importantly, 


25 


substantial  reductions  in  diffuser  losses  are  expected.  The  losses  are  discussed 
in  the  following  sections,  which  consider  two  all-supersonic  configurations.  In 
the  first,  a  single  large  diffuser  is  used  to  recompress  the  flow  prior  to  expan¬ 
sion  in  a  subsequent  turning  nozzle  array.  In  the  second  case,  a  diffuser  is 
machined  as  an  upstream  integral  part  of  each  noetic  element.  In  this  idealized 
ease  it  is  no  longer  necessary  to  turn  the  flow  i„  the  supersonic  expansion,  and 

a  laser  cavity  would  follow  each  array  of  diffuser-nozzles  fitted  to  a  straight 
duct. 


Gasklnetic  Losses  In  Nozzle  Diffuser  Systems 
The  nonequilibrium  losses  i  ave  been  calculated  for  representative  nozzle 
diffuser  combinations  using  both  a  two-temperature  and  a  four- temperature 
(Appendix  C)  model  for  the  N,,  -  CO,,  system.  Figure  8  is  an  example  of  the  latter 
procedure  applied  to  the  centerline  of  a  simple  symmetric  supersonic  configura¬ 
tion.  The  flow  expands  from  M  -  1.3  through  a  15°  half-angle  wedge  nozzle  to  a 
lasing  channel  which  has  ten  times  the  cross-sectional  area  of  the  nozzle  throat 
and  a  length  of  over  200  throat  heights,  A  7°  wedge  diffuser  then  recompress.s 
the  flow  back  to  equilibrium  at  M  -  1.3  prior  to  subsequent  re-expansion  and 
laser  energy  extraction.  The  figure  is  for  a  40%  N,,,  10%  CO,,,  and  50%  He  mixture 
at  To  -  2000%  with  pQh  -  1  atm  cm.  Tr  denotes  the  N,,  vibrational  temperature, 

Tv  and  Tl  the  upper  and  lower  vibrational  levels  of  CO,,  respectively,  and  T  the 
translational-rotational  temperature.  All  four  temperatures  are  seen  to  approach 
each  other  in  the  diffuser,  however  the  lag  in  ^  and  T^  is  such  that  an  extended 
throat  region  (1-2  throat  heights  long)  will  generally  be  required.  The  specific 
entropy  production  rate  (ds/dx)  may  be  found  using  the  relation 

dS  -  d!5  +  !i+da 

rjl  *  ip  T  Ijl  ^ 
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rfbe.,8  dQ  deN  -t*^eu  1  dc'L  +  ^LASE^ '  Here  dejj’  deu’  dnd  deL  reP resent  the  in- 

BR  ^  '■  :  !  ‘  '■  'l  .  J|  1 

crease  in  tie  various  vibrational  energy  storages  per  unit  mass  of  mixture,  with 

dQ  being  the  change  in  translational  rotational  energy  and  de^^,  the  energy  re¬ 
moved  oy  lasing.  The  lower  part  of  Figure  8  shows- that  the  entrbpy  foi^  this 
particular  flow  was  about  equally  produced  in  the  initial  nonequilibrium  expan- 
sion  and  in  the  lasing  channel,  with  a'  smaller  contribution  from  the  diffuser. 
Neglecting  changes  in  Tq,  the  ratio  of  pQ  between  exit  anfl  entrance  equilibrium 
states  ie  found  from  the  total  entropy  change  to  be  **  0.80. 

As  discussed  in  Section  III,  use  of  a  carefully  Shaped  combination  nozzle 

8  «  * 
and  lasing  channel  can  reduce  losses  of  the  total  pressure 4  .  Larger  scale  diffu¬ 
sers  would  accomplish  equilibration  near  the.  entrance  region  with  a  further 

r 

reduction  in  loss.  While  it  is  not  geometrically  possible  to  use  such  diffusers 

1  ,  r 

with  symmetric  supersonic  nozzle  elements  (such  as  Figure  7),  their  use  with 

«  ■ 

supersonic  turning  nozzled  could  rqduce  the  diffuser  Pq  loss  dhe  to  flow  nonequi¬ 
librium  to  1-2%.  The  calculations  show  that  with  advanced  nozzle  contbyrs,  , large 
diffusers,  and  optimized  running  conditions,  it  should  be  possible  to  build  super¬ 
sonic  configurations  with  less  than  10%  nonequilibrium  Pq  loss  without  significant 

change  in  X.  i  1 

1  s 

Fluid  Mechanical  Losses  i 

*  i 

Fluid  mechanical  loss  calculations  have  been  carried  out  for  a  represents- 

'  i 

live  series  of  diffuser-nozzle  combinations.  Air  was  used  £<pr  the  working  fluid, 
with  T^  2000  K,  PQh  -  1  atm  cm,  and  h  =  0-.  1  cm.  The  normal  cavity  Mach  number 
was  4.  This  flow  is  fed  into  a  2-D  diffuser  (area  ratio  =  10,  half-angle  =  7°), 
followed  by  a  0.1  cm  long  throat  and  a  15p  half-angle  nozzle.  t' 

,  >  i 

The  system  of  oblique  shocks  in  the  wedge  diffuser  produced  4  calculated 

_  'i  i  1 

Po  loss  of  ~  10%.  The  boundary  layer  was  assumed  to,  undergo  transition  at  the 


*  1 
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first  shock  boundary  layer  interacLi-jn.  An  integral  analysis  was  used  up  to  this 
point;,  followed  by  the  turbulent  boundary  layer  procedure  of  Sasman  and  Cresci. 

The  latter  was  initiated  with  the  assumption  that  the  momentum  thickness  (9)  was 
preserved  across  the  interaction.  Values  of  the  dimensionless  displacement  thick- 
ness  (6  /h),  equal  to  ~  1.0  and  0/h  ~  0.5  were  obtained  at  the  nozzle  exit  af^er 
carrying  the  calculation  through  the  individual  diffuser-nozzle  channel.  When 

i  ' 

mixed-  at  constant  area,  this  produced  a  pQ  loss  of  48%. 

Much  smaller  losses  were  found  when  the  diffuser  scale  was  increased  by  two 
orders-of-magnitude.  In  this  configuration,  a  single  large  diffuser  would  feed 
an  array  of  turning  nozzles.  The  diffuser  exit  boundary  layer  accounted  for  only 
~  of  the.  mass  flux  and  could  be  bled  off,  directly  reducing  X.  The  individual 
nozzle  boundary  layers  now  start  and  develop  in  a  favorable  pressure  gradient, 

,nnd  both  'analysis  and  experiment  show  the,.,  to  be  laminar.  Exlt  thicknesses  are  now 
found  to  be  <5*/h  s  0.8  and  0/h  «■  0.03,  resulting  in  a  mixing  loss  of  15%  in 

i 

p  . 

o 

Including  a  representative  nonequilibrium  loss  of  10%  in  p 

o’ 

the  collectfed  effects  lead  to  a  staging  figure  of  merit  X/ln(p  /p  )  of 

oi  oe 

40J/gm  for  ;the  turning  nozzle  configuration.  This  is  over  twice  that  for  the 
combination  diffuser-nozzle  elements,  and  is  considerably  larger  than  that  obtain¬ 
able  with  conventional  supersonic-to-subscnic  diffusers  and  nozfcles.  Further,  it 
is  probable  that  the  boundary  layer  will  relaminarize  in  the  small  diffuser  case 
due  to  the  high  acceleration  in  the  throat  region/18^  This  has  been  indirectly 
observed  in  laboratory  experiments,  and  it  should  make  the  losses  quite  similar 
for  the  two  supersonic  cases  studied.  Boundary  layer  control  through  bleed, 
mass  injection,  and  wall  cooling  are  being  studied.  The  diffuser  roust  have  a 
turbulent  boundary  layer  in  order  to  remain  attached;  however,  it  should  be 


pnillblti  to 


di'une  the  shock  losses  with  use  of  contoured  walls,  at  least  in  the 


large  scale  case.  The  practical  problem  of  diffuser  starting'  '  ifi  also  being 

-1  m  Bn  f'l  si >  *  ! ,  •  FSJ!  ■  i  ss 

considered. 
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TEMPERATURE  (°KxlO"3) 


Mgure  It  Vibrational  energy/mole  in  N2  normalized  by  the 
characteristic  energy  of  the  first  level  R6  (0  =  3350°K) . 
Ev  represents  the  total  vibrational  energy  while  Ev  *»  1 
represents  only  the  first  vibrational  level. 


mm 


7  inch 

ttwtV  T ub  c 


Figure  2.  Sketch  of  the  laser  optical  system  and  test  section 
used  to  measure  gain.  (i  =  diffuser  plate,  f  =  10. 6y  filter, 
bs  =  beam  splitter,  in  =*  mirror,  w  =  germanium  window). 


Side  view  of 
Nozzle  Slocks 


Detector 


Td 


Au-Ge 

Detector 


CO^  Gas  Laser 


(a)  Short  time  trace  oscillograph.  Vertical  Scales: 
Gain  (lOmv/division) ,  pressure  (250psi/division) ; 
Horizontal  Scale:  Time  (ltnsec/division) . 


Figure  3.  Typical  oscillographs  for  an  experimental  run.  Results  for  the 

5.18%  C0_  gas  mixture  are  shown.  Measured  gain  =  .55/moter  at 

T  -  2406°K  and  p  =  39.8  atm. 
o  o 
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A  Calculation  of  Mixed  Mode  Contributions 

to  Absorption  in  C02  at  10. 6ll 

•  1  .  I  j 

S.  A.  Munjee  and  Walter  H.  Christiansen* 

v  'WV'U  to''1  "■  ‘  •.  t 

Aerospace  Research  Laboratory,  University  of  Washington 
Seattle,  Washington  98195 

A  rise  of  temperature  in  a  molecular  system  causes  the  higher  energy 
levels  to  become  more  populated.  For  the  case  of  C02  absorption  (or,  in  a 

lasing  medium,  gain)  calculations  based  only  on  single  line  considerations 

*  < 

are  sometimes  inappropriate,  as  the  contributions  from  some  higher  energy 
levels  can  be  appreciable.  This  note  summarizes  this  effect  for  the  10. 4u 
band  of  CO^  The  calculations  are  compared  with  the  experimental  measure-  • 
ments  of  others. 

.  The  equilibrium  absorption  antributions  to  the  P-branch  00°1  t  10*0 (P2) 

transitions  due  to  near  resonant^ transitions  from  mixed  modes  of  the  type 
%  9,  * 

OX.  1  ^  IX  D  have  been  estimated,  where  X  is  any  positive  integer  and  i,  the 
angular  women turn  quantum  number,  has  appropriate  values.  For  the  range  of 
296°K  tc;  800°K,  it  was  determined  that  only  the  R-branch  transitions  01xl 
t  ll10(R7r)  and  02°1  Z-  12°0(RE)  contribute  significantly  to  the  overall  absorp¬ 
tion  of  the  P-branch  lines  because  near  resonances  along  with  appreciable 
population  of  energy  levels  exist  only  for  these  transitions.  For  example, 
the  02  1  ^  12  0(RA)  transitions  cannot  contribute  because  their  wavelengths 
are  not  sufficiently  close  to  those  of  the  PE  lines;  on  the  other  hand,  the 
transition  of  the  type,  002  £  lOl,  may  have  close  resonances,  but  the  popula¬ 
tion  of  those  states  are  too  small. 


*Assoc.  Prof.,  Dept,  of  Aeronautics  and  Astronautics 
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For  th<  calculations,  vibrationally  non— degenerate  energy  levels  of 

■ 1  ■  I  ■  EH 

F-f-  w  Hr 


,  *  *  i-  ^  f  1  ^  .  i>  (  •  p,'^'  \  ~  .  M  *  /. '  p 

molecules  of  a  gas  at  moderate  pressure  (collision  hroadened)  are  assumed, 


for  which  the  absorption  coefficient  at  line  center  for  a  single  spectral 


line  can  be  written  as:* 


VI, 


■  :  !«,:  K  ■  I 


. 1.^UL  e 

a  (A  )  «  -  i _ °fN  _  ISm  ‘tL. 

V  1  "  &■% 


a) 


where  the  subscripts  U  and  L  refer  respectively  to  the  upper  and  lower  levels 


of  the  rotational-vibrational  transition,  A^  is  the  wavelength  of  the  spec¬ 
tral  line  center,  is  the  radiative  lifetime  for  the  line,  the  g’s  ard  the 
rotational  statistical  weights,  and  the  N’s  are  the  population  densities  of 


the  rotational-vibrational  quantum  levels.  The  optical  broadening  collision 
frequency  Vc  is  Nacr,  here,  N  is  the  number  density,  a  is  the  mean  relative 
speed  of  the  molecules,  and  C  is  the  optical  broadening  cross  section. 

If  a  Lorentzian  profile  is  assumed  for  a  single  spectral  line  of  wave¬ 
length  A  the  absorption  coefficient  in  the  neighborhood  of  the  line  center 
may  be  expressed  as:2 


o(A) 


°<XUL  > 
o 


1  +  cfr<i  -  T^2 

VC  A  .  aul 


(2) 


Tha  total  absorption  coefficient  is  then  obtained  by  repeatedly  applying  Eq. 
(2)  to  contribu  ng  transitions  and  summing;  i.e., 

a(A}  =  apSCA)  +  a^CA)  +  2aR7r(A)  (3) 
Each  of  the  terms  on  the  right  hand  side,  of  Eq.  (3)  is  a  sum  of  contributions 
from  the  vibrational-rotational  lines  within  each  transition,  and  the  coeffi¬ 


cient  2  in  its  rightmost  member  accounts  for  the  two-fold  vibrational  degen¬ 
eracy  of  the  Tf  states.  [When  Eq.  (1)  is  applied  to  the  Rtt  transitions,  the  7f 
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States  dre  treated  as  though. they  were  vibrationpily  non-degenerate,  with 
their  populations  giver-  uy  the  Boltzmann  distribution.] 

In  order  to  evaluate  Eq.  (3)  [in  terms  of  Eq.  (1)],  the  radiative  life- 
time  for  spontaneous  emission  must  be  known.  This  -is  related  to  the  Einstein, 
coefficient  for  stimulated  emission,  and  can  be  expressed  ![Eq.  (2.17),  Ref.  3] 

in  terms  of  the  J  dependent  and  vibrational  (J  independent)  contributions  to 
the  transition  moment.  Thus,  ' 

SL 


tul  3h 


647^  .3SfJU)KUL 

(Aul  }  ~r~i~r 

o  BU 


(4) 


where  h  is  Planck's  constant;  S^),  the  rotational  contribution  is  taken 
from  P.  127  and  p.  208  of  Ref.  4,  and  is  the  matrix  element  of  the  dipole 
moment  squared.  To  obtain  the  radiative  lifetime  for  the  P2,  R2,  and  Rtt 
transitions,  is  treated  as  a  constant.  Such  an  approximation  seems  rea¬ 
sonable  because  these  radiaqive  transitions  do  not  involve  a' change  in  !Lhe  1 
quantum  number  of  the  bending  mode  of  C02.  The  numerical  value  of  (1.37 
x  10  39  erg-cm3)  is  obtained  by  applying1 Eq.  (4)  to  the  P(20)  line  of  the  PZ 
transition  with  tul  -  5.38  sec.5  Different  authors6'7  have  reported  a  weak 
variation  in  line  width  (hence  a  )  with  the  rotational  quantum’ number  ‘for 
the  PI  transition.  However,  variations  between  authors  show  wider  differences 
than  the  variation  with  rotational  quantum  number.  'Furthermore,  this  depen-  ’ 
dence  for  transitions  other  than  the  PZ  has  not  been  established.'  In  the  : 
calculations  therefore,  a  has  been  considered  a  constant  with  its  value  (1,3 
x  10  14  cm2)  taken  from  Ref.  8.  • '  '  ,  ! 

Finally,  the  various  wavelengths , (1^  )  needed  to  complete  the  evalua¬ 
tion  of  E,.  (3)  can  be  obtained  from  Eqt.  (5.1)  -  (5.4)  once,  the  rotational- 
vibrational  constants  of  the  molecule  are  known.9  ' 

PZ:  (V  >  1  =  Ec  +  ‘  yj(J+D]  '  (5.1; 
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I 


o 

Prr: 

i 

<V.  ,_1  ■  E„  +  SulJO-l)  -  i2J 

0  u 

(5.2^ 

_  i 

-  El[J(J+1)  -  +  (-1)J  lq  .(j-D  . 

-  2  u 

-  qLJ(J+l)J 

I 

RE: 

(XUL  ^  1  =  Eo  +  But(J+l)(J+2)]  -  Bl[J(J+1)] 

(5.3) 

'  ,  f ,  ‘ 

Ptt: 

(AUL>  1  =  Eo  +  V<J+l)(J+2)  -  Z}]  -  B  [J(J+D 

O  Aj 

-A*] 

'*  ililifqyCJ+l)^)  -  qLJ(J+l)]  ( 

i.  i  '  _ 

•Ini  Eqs*  (5-1>“  EQ  is  the  difference  between  the  vibrational  energy 

!  ‘  !  ^EvibV°f  the  upper  and  lower  levels,  expressed  in  cm"1,  the  B's  are  the 

rotational  constants,  and  the  q*s  are  the  A-type  doubling  constants,  both  of 
whibh  are  expressed  in  cm"1;  J  is  the  rotational  quantum  number  of  the  lower 

V  ♦  V,  • 

level.  The  calculation  of  the  B's  and  q's  cannot  be  readily  accomplished 
because  the  lower  levels  of  each  of  the  3  transitions  are  in  Fermi  Resonance 
with  other  levels  of  the  molecule.  For  the  PE  transition,  the  constants 
'Q  ’  WerC  taken  'fro!n  Ref  •  10=  the  RE  transition  from  Ref,  II,  and  for  the  R* 

transition  from  Ref.  12.  This  is  summarized  in  Table  I.  Using  the  values  in 
able  I  and  Ejs.  (5.1),  (5.3)  and  (5.4),  "near  resonances"  were  calculated; 
these  are  shown  in  Table  II  for  some  selected  JL  values  of  the  PE  transition. 

I  Calculations  of  the  total  absorption  coefficient  [Eq.  (3)]  for  the 
wavelengths  of  the  PE  transitions  were  carried  out  for  values  of  J  from  14 

i  L 

,  through  28.  Because  of  interest  in  the  strong  P(20)  laser  line,  results  for 

i  thiS  Wavelength  are  shcwn  ^  Fig.  1.  This  figure  compares  the  results  of 
calculations  made  in  this  study  with  the  experimental  measurements  of  Refs. 

14  and  15.  The  broken  curve  represents  c^,  the  single  line  center  absorption 
coefficient  for  the  P(20)  line  only.  The  solid  curves  represent  ct,  the 
t0tal  absorPti°n  coefficient,  which  is  comprised  of  aQ  and  contributions 
from  mixed  modes.  Interest  has  persisted  in  the  Rtt(23)  line  and  the  PE(20) 
line  as  these  two  are  believed  to  be  very  nearly  resonant.  The  authors  of 
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a  wave- 


Ref.  14  believe  this  difference  to  be  within  several  hundredths  of 
number;  in  this  study  (Table  II)  the  difference  is  estimated  as  .011  cm"1. 

In  the  lower  temperature  range  (296°K  -  425°K'  the  calculation  of  a  for 
p  »  1/3  atm.  is  compared  with  the' experimental  measurements  of  Ely  &  McCubbin.14 
While  there  is  considerable  scatter  in  the  data  Figure  1  indicates  that  the 
calculated  values  of  a  are  generally  somewhat  lower  than  those  measured. 

.  However,  if  the  uncertainty  of  ±8%  indicated  in  Ref.  14  is  assumed  to  apply 
throughout  the  lower  temperature  range,  the  agreement  is  better.  The 
contributions  from  the  mixed  modes  at  300°K  is  seen  to  be  about  6.4%.  Ref.  14 
estimates  it  as  3.2%  based  on  the  assumption  that  the  line  intensity  ratios 
depend  only  upon  the  Boltzwac.u  factors  and  the  (2J  +1)  degeneracies.  If 
allowance  is  made  for  the  two-fold  vibrational  degeneracy  of  the  7r  states, 
the  estimate  of  3.2%  will  be  doubled.  Also,  Eq.  (2)  of  Ref.  14  does  not  seem 
to  allow  for  the  temperature  dependence  of  the  rotational  populations. 

In  the  higher  temperature  range  (above  400°)  the  calculations  of  a 
(p  =  1  atm)  are  compared  with  the  measurements  of  Leonard.15  Although  Leo¬ 
nard's  data  exhibit  a  fair  amount  of  scatter  the  figure  indicates  that  the 
total  absorption  (a)  is  considerably  greater  than  the  single  line  center  ab- 
sorptJ  jn  (ao).  At  400°K,  for  instance,  the  calculated  contribution  is  about 
17.7%  and  at  620°K  it  is  42.1%.  The  figure  also  suggests  that  the  measure- 
aents  correspond  more  closely  to  a  than  to  aQ.  Calculations  of  a  for  PE(14), 
PZ(16)  and  PE(18)  have  been  compared  with  Leonard's  data15  and  it  was  concluded 
that  the  calculations  agreed  with  his  data  to  within  10%.  Without  the  mixed 
mode  contribution,  the  difference  was  as  much  as  four  times  that. 

Figure  ?  shows  the  qt.antity  Si  «  (a/ao  -  1)  a3  a  function  of  pressure 
(for  seme  selected  temperatures)  for  the  case  of  PE (20).  The  results  for 
other  lines  are  similar  in  trend;  SI  tends  to  increase  with  pressure,  then 
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levels  off,  and  Increases  once  again.  .For  very  close  resonances  such  as 
the  PS (28)  -  Rir (12) ,  C^)T  "  0  for  pressures  approximately  .05  atm  to  1  atm. 
Wii-s  characteristic  is  shown  in  Figure  2  as  the  two  Broken  curves  (for  T 
“  296°K  and  T  *  800°K) .  For 'wavelength  differences  greater  than  5%B 
(B  s  .39  cm  l)  there  is  no  extended  region' of  pressure  for  which  this 
characteristic  is  observed. 

The  absolute  magnitude  of  the  mixed  mode  contributions  are  questionable 
primarily  because  the  rotational-vibrational  constants  of  the  molecule  are 


not  known  with  adequate  precision.  Despite  this,  the  behavior  of  ft  can  be 

expected  to  be  similar  to  that  in  Figure  2.  This  is  so  because  for  a  given 

value  of  "near  resonance"  [1  «  ~  ~  in  Eq.  (2)]  the  effect  of  the  un- 

'  AuLo 

certainty  in  A  can  be  offset  by  .  compensating  change  in  V„,  i.e.,  a  dif- 

w 

ferent  pressure.  This  suggests  that  a  sensitive  experiment  conducted  to  i- 


dcntif'jr  would  enable  setter  values  of  the  Xotational-vib rational  con¬ 

stants  to  be  inferred. 

In  review,  it  is  believed  that  the  inclusion  of  the  contributions  of  the 


mixed  modes  gives  satisfactory  agreement  with  experiments  on  absorption  at 
10.6]l,  notwithstanding  the  assumptions  regarding  a  and  K^. 

The  authors  are  gratei^l  to  Dr.  G,  A.  Tsongas  for  his  very  helpful  com¬ 
ments  and  suggestions.  This  work  was  supported  by  the  U.S.  Army  Research 
Office  (Durham). 
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CAPTIONS  AND  NOTES  FOR  FIGURES 


Figure  1:  TOTAL  ABSORPTION  AND  SINGLE  LINE  CENTER  ABSORPTION  FOR  PE (20) 

V  / 


X  1/3  atm.,  Ref.  14.;  61  atm. , 0  1/2  atm.,  from  Ref.  15. 


The  broken  curve  is  cy.  The  two  solid  curves  are  calculations  of 


a;  the  lower  temperature  range  (296°K  -  425°K)  is  a  calculation 


for  p  =1/3  atm. ,  the  higher  temperature  range  is  for  p  -  1  atm. 


Figure  2:  CONTRIBUTION  FROM  OTHER  LINES  TO  SINGLE  LINE  CENTER  ABSORPTION 


. 

The  solid  curves  are  for  PE (20)  with  a  =  1.726xl0"3  cm"1 (296°K) , 


1.243x10  2  cm  J(500°K),  1.766xl0"2  cm-1 (300°K) .  The  broken  ones 


are  for  PE(28)  with  c*o  -  I.l4xl0"3  cm"1(296°K)  and  1.858xlO"2cm"J 


(800°K) . 


-  *3 
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COMMUTATION  OF  POPULATION  INVERSION  IN  C0o-N„  USER 

2  2 

MUIilRii  Brit  I  *  t  FIKT1*  TUMEfl  t  „Em« 
Oktay  Yesil 


INTRODUCTION 


The  invention  of  the  laser  is  perhaps  one  of  the  most  important  tecnnical 
developments  of  the  last  decade.  This  device  has  given  scientists,  and 
engineers  a  new  tool  for  generating  coherent  radiation.  Their  present  use 
in  many  areas  and  their  future  applications  have  been  well  recorded.1’2 

There  has  been  much  theoretical  work  done  on  the  kinetic  processes  of  relaxa- 
tional  phenomena.  In  some  of  their  previous  works,  Rubin  and  Shuler3 > 4 > s,  and 
Montroll  and  Shuler6  have  made  a  theoretical  study  of  the  collisions!  and  radiative 


relaxation  of  a  system  of  harmonic  oscillators  with  an  infinite  number  of 
levels.  Using  the  Boltzmann's  formula  and  under  the  assumption  that  the 
quantum  energy  is  very  much  less  than  kinetic  energy  (i.ts.,  hv/kT  «  1), 
they  obtained  the  vibrational  relaxation  master  equations.  They  found  that 
the  detailed  relaxation  behavior  depends  strongly  on  the  form  of  the  colli- 
sional  probabilities  and  that  radiative  transitions  do  not  change  the  relax¬ 
ation  behavior  of  the  system. 

Basov  et  al  7,  on  the  other  hand,  have  carried  out  the  kinetics  of 
systems  by  mode  analysis  rather  than  on  an  individual  level  basis.  Using 
mode  analysis,  many  authors  have  analyzed  the  population  inversions  of  CO^- 
^2  laser  mixtures  in  rapid  expansion  nozzles  using  the  generalized  energy- 
level  diagram  (fig.  1)  and  have  calculated  the  gain  between  the  levels  (00°1) 
(10°0)  of  C02,  1  =  10.6  micron. 

To  the  best  of  our  knowledge,  the  possibility  of  obtaining  gain  between 
the  levels  (02°0)- (0110)  of  CO^,  16. by  (simultaneously  when  lasing  the  CO^ 
molecule  at  10. 6y)  has  not  bean  reported  so  far.  Hence,  this  work  is  under¬ 


taken. 


SYMMETRY  BENDjNG  ASVMMETftEC 

VIBRATION  C  **1  VIBRATION  r 


FIGURE  1.  GENERALIZED  PARTIAL  ENERGY  DIAGRAM  OF  C0o-N-  MIXTURE 


This  work  is  devoted  to  the  establishment  ot  the  conditions  umier  Ihich 
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sutrtcient  gain  between  the  levels  CC>2  (OOpl)-CO  <10e0),  10.6yr  as  well  as 
between  CO^OVOycO^O),  16.18b,  is  obtained  In  CO^-He  laser:***  fixtures, 
varying  gas  concentrations  and  the  -atios  of  vibrational  to  translational  tem¬ 
peratures.  Initial  number  densities  "of  all  vibrational  levels  were  taken  as 
rero  except  CO^OO’l)  and  Njfv-l),  wt  ch  were  given  a  Maxwell-Boltsmann  tem¬ 
perature  of  3000-K.  The  gas  was  held  at  a  constant  pressure  of  1  atmosphere 
U4.7  psi).  Because  of  the  complexity,  a  finite  number  of  levels  of  harmonic 
oscillators  has  been  considered,  and  the  vibration  to  vibration  (V-V)  transitions 
have  been  neglected.  Actually,  since  V-V  transitions  are  faster  than  the  others, 
we  have  investigated  the  effects  of  their  absence.  Practically,  this  can  be 
achieved  by  reducing  ^-concentrations.  Therefore,  we  dealt  with  low  COj- 

concentrations  only  (up  to  10*  by  molar  fraction)  in  the  calculations. 

THE.C02-N2  system 

C02  is  a  linear  triatomic  molecule;  It  has  four  vibrational  degrees 
of  freedom:  the  symmetric  stretch  mode  ^  the  doubly-degenerate  bending 
“°de  V2'  and  th«  asymmetric  stretch  mode  vy  Nj  is  a  symmetric  diatomic 
molecule.  It  has  only  one  vibrational  degree  of  freedom  which  is  commonly  ' 
designated  by  v.  (see  figure  1). 

Be  have  considered  a  CO^-Nj  system  for  which  the  kinetics  are  well  - 
known"- 5  and  which  was  initially  excited  to  a  nonequilibrium  vibrational 
distribution  by  different  methods"-7-".  In  our  work,  the  levels  of  the 
system  have  been  allowed  to  reach  a  quasl-steady-state  condition;  then  strong 
lasing  at  10.6b  has  been  turned  on  by  q-switching  a  laser  cavity.  This 
quickly  bleaches  the  CO^OO7!)  to  COjdO’O)  transition,  which  increases  the 
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C02<02  P°Pulation,  through  Feml-resonance,  to  give  a  population  1 

inversion  between  CO^OZ'O)  pod  CcyoUo)  levels.  Here.we  have  used  the 
level  analysis  rather  than  the  mode  analysis.  It  cannot  be  used  since  we  hope 
to  get  an  inversion  within  a  mode  (i.e.,  negative  temperature)  during  lasing,  and 
therefore  the  mode  is  not  in  equilibrium  as  the  mode  analysis  assumes.  ’ 

modet:  _  ] 

Figure  2  shows  the  newly t developed  schematic  model  of  the  major 
participating  vibrational  energy  levels  for  CO^  laser  gas  mixtures. 

This  was  done  by  writing  down  all-  the  possible  levels  in  the  order  of  '  ' 

i  ' 

■  increasing  value  of  their  vibrational  energies  until  the,  cut-off  line 
(strong  dashed  line  at  E-3775  cm'*)'.  Energy  is  measured  from  the  ground  state  of. 
C02(00-0)  and  SjM)  and  only  the  pertinent  vibrational  levels  (totalling  25)  have 
been  shown.  The  effect  of  levels  above  the  cut-off  line  has  been  .neglected  because 
the  higher  levels  have  many  fewer  particles  af  this  vibrational  temperature. 

The  population  of  levels  becomes  redistributed  as  the  temperature 
changes  according  to  the  Boltzmann’s  formula:  , 

I 

n.  .  ,  *  , 

1  +  1  ,  4E,  .1  1  & 

n  pIT  1  (1) 

i  v  i 

where  AE  is  the  energy  difference  between  the  levels  i'+  1  and  i  of  a  single  1 

mode,  i  .  . 

Each  of  the  vibrational  levels  shown  is  split  into  a  series  of  sublevels 
due  to  the  quantized  rotational  energies  of  the  molecules  as  a  whole,  but 

f  '  !  ' 

these  are  not  indicated.  Table  I  shows  the  vibrational  energy  levels  of  1 
C02-N2  mixture.  Column  1  gives  the  levels  and  Column  2  shows  the 
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TIGURE  2.  SCHEMATIC  MODEL  FOR  A  (CO.-NJ  SYSTEM 
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corresponding,  energies  in  cnT1,  Column  3  gives  the  multiplicity  (or  degree 

of  degeneracy)  of  the  levels.  This  multiplicity  occurs  only  for  v10. 

n  the  model,  {  }  s  in  fig.  2  contain  the  levels  which  have  nearly 

eC*U  energies  and  each  {  }  has  been  labeled  as  B  (Box  ),  B  ,  B  ,  etc. 

1  JL  2 

Hence,  fast  transitions  appear  to  justify  a  vibrational  model  which  groups 
.  the  participating  levels  into  11  boxes.  It  is  assumed  there  Is  local 
equilibrium  within  each  box  but  not  between  boxes. 

In  this  work,  only  the  following  relaxation  processes  have  been 
considered1 1 : 

1.  -Translation  to  vibration  (T-V),  strong  solid  lines 

C02(v1v2v3)  +  M  t  C02(vlv2-iv3)  +  M 

2.  Special  transitions  (SP) 
a)  SPx,  light  dashed  lines 

C02(V1V2V3)  +  N2(V)  Z  C02(V1V2V3  +  2)  + 
k)  Usht  long  dashed  lines  • 

C°2(V1V2V3)  +  M  *  C02(Vi  V3  V1}  +  M 

3.  Radiation  transitions  (RAD),  wavy  lines 

Here,  M  Is  called  the  collision  partner  and  M  =  C02,  N2  or  He. 

We  note  that  (V-V)  transitions  have  been  neglected.  Because  the 

probability  of  them,  ^V-V  ~  irf’3^12^  i  ^  . 

*  ^  i«e.f  they  are  about  10  times  faster 

compared  to  the  other  transitions  in  our  work. 

GOVERNING  EQUATIONS 

1.  Master  vibrational  relaxation  equations 

formally,  any  simple,  linear  relaxation  process  can  be  described 
by  an  equatior  * 


u 


u 


u 


5?  '  -  k  <U-De> 

Where  “  13  ,Ua"tity  ”hiCh  ls  -d  »•  is  its  equilibrium  value. 

T  IS  the  charactaristic  tine  of  the  relaxation  process. 

The  naster  vibrational  relaxation  equation  for  the  nurber  density, 

■>t  at  level  i,  contains  the  tarns  of  all  possible  relaxation  processes, 
i.e.. 


dt  ;T 


dni,  dni  dn 

dt  T-V  +  dt  \-V  +  dt — \-V 


dni 

+  —) 
dt  'RAD 


(intra  mode)  (inter  mode) 

Neglecting  (V-V)  transitions  and  including  the  special  transitions 
and  introducing  the  fraction  of  the  level  i.  ^  ,  ls  the  total 
number  of  particles)  we  have: 


dx±(t) 

"IT 


dx, 


dx, 


\  _  i j  dx 

T  dt  T-V  +  “dt^SP  +  dt^SAD 


For  the  complete  set  of  master  vibrational  relaxation  eouations  used  in 
the  computer  program  see  the  Addendum. 

2.  Equations  for  the  Gain  Calculations 

From  reference  (13),  the  expression  of  the  line  center  gain 
coefficient  for  a  collision  broadened  line,  is  given  by 


\  2 

a  =  — ■ — A 
4ttT(t:Av  ) 


2  nj 
gj 


V 

gj" 


•) 


where  to  is  the  line  width  and  is  related  to  the  collisional  frequency 

V  by  V  =  7rAv  . 
u  c  c 

The  expression  for  V  ,  in  the  case  rn  „  „ 

c»  j-u  cne  case  CO^^-He  mixture,  is: 
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13, 


V  e  N  fcT  c  + 

C  wco2iaco2cco2  + 


He 


N„ 


N, 


CO, 


au  CIT  + 

He  He 


N, 


CO, 


\  V 


(3) 


where  N's  are  the  population  density  of  the  species,  a4 s  are  the  collision 


cross-section  of  each  species  with  C02> 


CO 


*  1.3  x  10“ 14 cm 2 


He 


N, 


1.7  x  10“ 1  'em 


8.7  x  10“15cm2 


C’s  are  the  mean  relative  velocity  of  the  species,  c  *  /§££• 

*  i  VTry 


y, 


— -  +  — 
t-i  Di- 


X  is  the  transition  wavelength,  and 

.  10. 6u  for  the  radiation  between  C0_(00°1)  -  CO  (10°0) 

X  *»  {  *  2 

16.18y  for  the  radiation  between  C02(02°0)  -  C02(0110) 

T  is  the  radiative  life  time  for  the  transition,  and  from  reference  (14) 

T  .  ^-4.17  sec.  for  the  radiation  at  10. 6y 
2.08  sec.  for  the  radiation  at  16.18y 

Ej  is  the  number  of  particles  in  the  Jth  rotational  band  of  the  upper  level. 
nj  *s  t*ie  number  particles  in  the  J'th  rotational  band  of  the  lower  level, 
gj  is  the  rotational  degeneracy,  and  g  =  2J  +  1 

J 


8J" 


2J"  +  1 


We  note  that  for  the  radiation  at  10. 6y,  P-branch  transitions  and,  for' 

the  radiation  at  16.18y,  Q-branch  transitions  give  the  most  gain. 1 0 ’ 1 3 * 1 4 

Furthermore,  _J  and  J  in  equation  (2)  have  been  evaluated  in  the 
8J  8J" 

following  way  for  the  case  at  16.18y: 
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_ _ •  n(02e0)  M  (upper  level) 

gj  gJn(02°0)  N  N 


_ •  •  n(011Q)  „  (lower  level) 

n  rm  ir\  w 


gj-*  gJ»n(0110 


On  the  other  hand,  from  ref.  (12), 


gjn(02°0) 


gJ.n(0TT0) 


_1_  exphJ(J  +  1) 

Q  i 

1 


exPf"  J‘f(J>  +  1)  |£] 


'R 


•  T 

where  QR  A  :  rotational  partition  function 


hefi 

gr  =  — =  0.56  :rotational  characteristic  temperature 
B  a  0.39  :  a  coefficient  dependent  on  the  vibrational  level 


J  -  J 


Max 


_[  kT 
i  2hcE 


corresponds  to  the  peak  gain 


(4) 


(5) 


Ke  note  that,  ^  ^  and  -  come  directly  from  our  computer  program. 

For  more  detail,  see  references  (10,  13,  and  14). 


RESULTS  AND  DISCUSSION 

To  obtain  a  solution  for  a  given  set  of  initial  conditions,  the  coupled- 
equations  must  be  numerically  integrated.  In  the  set  of  master  vib rational 
relaxation  equations  (Appendix  A),  the  independent  variable  is  time,  t(  or 
in  nondimens ional  form,  t  =  t  Ek^^g  ,  see  the  Addendum).  There  are  several 
dependent  variables:  the  rate  constants  (see  Appendix  A)  which  have  been 
obtained  from  experimental  results  given  in  reference  (16)  at  the  correspond¬ 
ing  translational  temperature,  and  the  total  number  of  particles,  which  has 
been  calculated  with  the  state  equation  (i.e.,  P  =  NkT) ,  for  the  corresponding 


61 


translational  temperature,  T.  TUn  the  populations  (or  the  population 
fractions)  of  the  levels  are  computed  versus  time.  The  10. 6y  small  signal 
gain  coeffi  :ients  have  bean  calculated  at  t'  ■  5  (time  to  reach  a  quasi- 
steady-state  condition),  and  t;ne  16.18y  gain  coefficients  at  t'  *  5.02 
(just  after  the  lasing).  The  effect  of  the  absorption  from  the  level 
C02(00°0)  to  C02(0.110)  on  the  gain  has  been  analyzed  and  found  to  be 
unimportant  for  the  C02(02°0)  to  CO^O^O)  transition  at  16.18y. 

The  system  of  equations  has  been  solved  in  FORTRAN  IV,  on  the  CDC 
6400  computer  (at  the  University  of  Washington)  using  a  5th  order'Runge- 
Kutta  technique,  subroutine  RKINIT*1 . 

TMs  technique  uses  error  estimates  based  on  the  comparison  of  tvo 

4th  order  solutions;  5th  order  results  are  obtained  by  an  extrapolation 
process. 

Some  typical  results  have  been  shown  in  figures  3  through  8,  and 
tables  II  and  III  for  the  following  data: 

P  *.  1  atm  :  total  pressure 

Tv  3000°K  •  initial  vibrational  temperature 

Vl  =  100  x  kf  x  N  :  for  the  laser  radiation  only  between  C02(00°1) 

and  C02(10°0) 

with  variable  Tv/X  ,  X^/X^  and  mixture, 

we  note  tha  :,  the  laser  radiation  can  be  produced  faster  than  any  other 
relaxation  process  by  a  higher  intensity  light.  Therefore,  we  have  assumed 

*1Boeing  Library  of  Mathematical  Routines 


it  is  100  times  faster  than  the  biggest  rate  constant  in  the  system. 

Figure  3  shows  the  variation  of  the  population  f  rac  w.ions  divided  by 
their  degeneracies  of  the  levels  with  the  noudimensional  time,  t""  for  the 
mixture  5,4C02  +  70%N2  +  25%He  and  Tv^  *  15,  The  population  fractions  of 
the  levels  C02(0110),  C02(02°0)  and  j?2  (v^O)  increase  while  the  levels  C02(00°0), 
C02(00°1^  and  N2<v=l)  eventually  decrease  with  t'  up  to  a  quasi-steady-state 
condition  (t  =5) .  At  t  =5.02,  however,  just  after  the  laying  has  been  turned 
on,  the  populations  (particularly  of  the  levels  C02(00°1)  and  C02(02°0))  are 
affected  strongly;  the  levels  C02(00°1)  and  C02(02°0)  are  bleached  immediately. 
There  becomes  the  possibility  of  obtaining  a  population  inversion  between  the 
levels  C02(02  0)  and  C02(0110)  because  their  relaxation  time  is  much  slower 

than  the  radiative  relaxation  time  of  C02(00°1)  and  C02(02°0)  transition. 

! 

These  curves  have  the  same  trends  and  almost  the  same  characteristics  for 
the  other  conditions. 

> 

ligure  4  shows  the  variation  of  the  population  fraction  over  degeneracy 
versus  the  energy  of  the  level  for  the  mixture  5%C02  +  70%N2  +  25%He  and 
*v/T  =  10  at  t  =5.  The  level  C02(00°1)  is  characterized  by  the  vibrational 
temperature  Tv  (solid  straight  line).  The  levels  C02(01l0),  C02(02°0), 
C02(0310),  and  C02(04°0),  however,  are  characterized  by  the  translational 
temperature,  T  (dashed  line).  Because  of  exclusion  of  (V-V)  transitions, 
this  line  has  a  slightly  broken  behavior  rather  than  being  straight.  We 
note  that  the  mixed  modes  (for  example  CO^IJ^O),  C02(021)  etc.)  have  be* 
spread  in  the  region  between  these  two  lines.  It  can  be  easily  seen  that 
we  can  never  get  population  inversion  between  the  levels  C02(02J0)  and 


(X^COl/O)  unless  with  the  lasing  action  at  10.6)1. 

Figure  5  shows  the  variation  of  the  gain  coefficients,  a  (both  at  10. 6y 
and  16.18y)  with  Tv/,j,  for  a  fixed  He- concentration  of  25%  and  with  the  ratio 
of  CO^  to  N2  as  a  variable.  Higher  C02  to  ratios  and  higher  vibrational 
to  translational  temperature  ratios  give  higher  gain.  At  Tv/T  **  7.5,  which 
corresponds  to  I  =  400°K,  no  16.18y  gain  has  been  predicted.  Calculations 
have  also  shown  that  the  gain  at  10. 6y  ir  higher  than  the  gain  at  16.18y  at 
these  same  cn.iditions.  The  curves  almost  have  straight  line  trends. 

Figure  6  shews  the  variation  of  a  with  X  /X.T  for  a  fixed  Fe-concentra 

C02 

tlon  (25%)  and  with  translational  temperature  as  a  parameter  (T  *  200,  300 
and  400°K).  This  figure  supports  the  conclusion  drawn  in  fig.  5.  It  is 
shown  that  higher  C02  concentrations  give  higher  radiation  (or  absorption) 
due  to  the  higher  density  of  C02  particles. 

Figure  7  shows  the  effects  of  the  mole  fraction  of  He  on  the  16.18y 
gain  at  fixed  XCQ  /X^  (=  1/14)  and  T  *  300° K.  Higher  He-concentrations 
(up  to  10%)  increase  the  gain  very  fast;  but  for  fractions  higher  than  10%, 
they  decrease  the  gain  slowly.  At  zero  He-concentration,  no  gain  has  been 
predicted. 


THE  LEVELS  /  DEGENERACY  ,  x„/g 


3000° K  ,  P  -  I  ATM  ,  Tv/T  =  15 
=  5%C02  +  70%  N2  +  25%  He 


Tiiwir  t 


FIGURE  3 


TABLE 


II 


u 


P  «  1  ATM. 
Ty  =  3000°K 


MIXTURE 

VT  . 

GAIN  COEFFICIENT,  a  (cT1) 

C02(02°0)-C02(0110) 

X  ■=  16.18m 

co2(oo°i)-co2 (10°0) 

X  -  10.6m 

5%C02 

15 

4.57 

5.21 

70%N2 

10 

0.957 

3.32 

25%He 

7.5 

-0.816 

2.35 

2.5%C02 

15 

2.396 

2.716 

72.5%N2 

10 

0.557 

1.697 

25%He 

7.5 

-0.413 

1.12 

i%co2 

15 

1.008 

1.1 

74%K2 

10 

0.242 

0.689 

25%He 

7.5 

-0.1603 

0.477 
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RATIO  OF  THE  INITIAL  VIBRATIONAL  Tv 
TO  TRANSLATIONAL  TEMPERATURE  ’  T 


FIGURE  5.  VARIATION  OF  OaaN  COEFFICIENT  WITH  T  / T 

v 
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RATIO  OF  CARBON  DIOXIDE  TO  X 
NITROGEN  MOLECULAR  FRACTION  ’  X 


FIGURE  6.  VARIATION  OF  a  AT  16.1 8u  WITH  X  /X 


MOLE  FRACTION 


MIXTURE 


!*  . 

Ua  */  itaI  itmo 

v«.  mw)  tm  u  j  « 

» 

6.65%C02  +  93.35%N2 

— i - - — ! - - — - - — ■ — 

0 

6.35%C02  +  88.66%N2  +  5%Ke 

.  5 

6%C02  +  84%N2  +  10%He 

10 

.5.665%C02  +  79.335%N2  +  15%He 

15 

5%C02  +  70%N2  +  25%He 

25 

4%d02  +  56%N2  +  40%He 

40 

2.5%C02  +  35%N2  +  62.5%He 

62.5 

MOLE  FRACTION  OF  HELIUM  ,  %  (by  volume) 

■a 

FIGURE  7.  EFFECT  OF  He  CONCENTRATION  ON  GAIN  AT  16.18y 


addendum 


Master  Vibrational  Relaxation  Equations; 


Consider  the  level  CO^Ol^O).  It  has  2  gains  and  2  losses  in  (T-V) 


transitions. 


Reactions:  0110  +  M  00°0  +  M 

H 


K 

020  +  M  -r  01*0  +  M 


M  »  C02(0110),  C02 (02” 0) ,  C02(10°0),  etc. 


Then 


» 

4»  •  •  •  • 


•  •  •  • 
+  •  •  - 


» 

» 

<-\j>  +  *.  w~Mw4w->  ** 

+  (-  kf,^.0Ni  'Wo)  +  ■.-»?*  W«t :M)~ltit_r*,*NNCr.D)) 
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u 


Or 


dt 


t;h,j  l  vW  iv.*  -Vf  I  “■) 

■  *  (■C/'f/'s  •‘os*!) 


Since 

X  N’ 
i*i  1 

=Nc°i 

Z  Nw.  =  Nfj 
><  NJ  N* 

and 

kh-H 

=2  kf 

+1-90 

lie  «  3  k r 

'  T1 

and 

,  etc.. 


V 


etc. 

from  the  following  relations 


/y 

{<JNv\1o)x 


Net fo 


Noo\ 


N0 i'o  JZ-Ncdb  .  c~6h  * 


1  d/T 

=  -L  e  1 


k<^-t  *  i  k< 


kUi~*i  =  (  i+i)k<^0 


u 


-  &£, 


kl-4-+i  1^0  e  ^  sC-t+Ok^Q  e  *f 


Now,  introducing  the  population  fractions  of  the  boxes  (B^  B2>  .  .  see  Fig.  2), 
i.e.,  X±  (=  -— ) ,  finally  we  have 

-777=  Z  kh-lM[(-*  eM-e/r)X,-X2)  +(J/2X3  -S  exp(-e/T)xJj 

Likewise,  doing  the  same  type  analysis  for  each  box  and  for  each  relaxation 
process  (i.e.,  T-V,  Radiations,  Special  transitions),  the  complete  set  of  master 
relaxations  equations,  used  at  the  computer,  are  given  as: 
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THE  SET  OF  MASTER  VIBRATIONAL  RELAXATION  EQUATIONS: 


J  U=  ^  H  [xi  '2  exp('  ®A)  XJ  *  hf  N  &  x, .  -x,  x10)+^  i^-xj 

i  ^  \  p 

71  ~JT~  ^  ^->£7NL^  eM-°A)  Xt-X2)+(3/2.  Xj  -3  exp(-£>/r)  X,)] 

+  ^  ^  (X7  X<1  “X*  XY0)  f  K2  Ig  (  f/4  Xj  -X^K,  I3 (%  XJ  -X2) 

+  K4  ^  “X2J 

~7r)r  Z  kf’-*°H  {(3  exp(-  ^T)  \~3/z  %  exp(-%)  Xjj 

+  kf N  [■34('xJx„-xJ  x,0|f(KtiifK,i<Xxtf..^xJ) 

+  K*T*  fXi  '  ^  Xj)  +(x2  -%  XjJ  Kj 

7?i=  1 0  [(%  exp  W  x3  -  %  X4J  *(*$&  Xfi  -  $ 

+  2  ki<r«0M  [x*  - %  fiXpC-®/rj  x4]  *  m  (x7  - «  x,) 

d7^  =  2  k-f,  ”,  M[  x7  -2  e*P(-e/T)  X^  +  kf  N  (x,  X,,  -  X ;  X„j 

*z\oZ»? 

+  (H^,i(X'4xrx5) 


dX  \ 

IUr  =  ^  k^v,  p/j  exp(-e/T)X4-2%x6J+(^/j  X*  -  H/3  «p(-ff)X 

'  [f7  eXp("®/<r)X6  +  K'-T<-  (3/4  X3  -  Vi  X6  ) 


~1 

6 


jr)T  "  2  ^  "  0M  ^ze*p(-e/-r)xs-x7)  +(3/2  x3  -3  exP(-  */r)x7)] 

+kf  w(x2  x(0  -  x7  x„|*  e*p(-%-)x6  -xJ 

+  Kt.Ii(V3Xjt-Xj) 


^  =  Z^'loM  |(^v(-«/Tjx7-%  Xj)j+^«  |%(xJX„-XJ  X,)] 

+  z  fcW«&  [%  "v  W  x»  -X3]  (fc  Xt  -%  Xj) 

TTxj-  =  (x„  ~Xu>)+kf  H  f(x^  +X7+%  xs)xti-(x,+X2+%Xj)Xw 

(^2)  ~ 

=  2  kft m  (x,a-x(,J  Vf(Mx2  +  %  X3j  X,„ -fe  xx7  +34  X3) X,] 
r  .  Ml)  .  J 


WHERE  Xj's  ARE  THE  POPULATION  FRACTIONS  OF  THE  BOXES  IN  FIG.  :! 


AT  THE  EQUILIBRIUM: 


Q 


where 


2  k,  M  M 
Ti-*0 


k-  C°2  >,  I .  Hz 


U->0  Nc02  +  Nn2  +  Uull  N“* 


O 


u 


2  k(<Z  m =kf<c°o  +  kfl"*  **  < 

Ml)  l/Vt»  CA/j.)  C^J 

and  they  are  related  to  (T-V)  transitions. 

i  co2 

q  rate  constant  that  a  C02_molecule  in  quantum  level  (01*0) 

wil.l  suffer  a  collision  with  another  CC^^olecule  and  make,  a 


transition  to  ground  level  (00° 0)  per  second. 

V 


r.  Afe 

fl-yO  ra^e  constant  of  the  collision  a  CO,,— molecule  with  a 


N  -molecule . 

kC°2 

■fl-+0  tbe  ra^e  constant  that  a  N_-molecule  in  quantum  level  (v»l) 

(>J2)  2 

will  suffer  a  collision  with  a  CO^-molecule  and  make  a  transition 

to  ground  level  (v=0)  per  second. 

Likewise  for  the  others. 

NC0  >  V’  NHe  are  the  total  nujnber  particles  of  CO  ,  N  ,  He,  respectively. 
£  £  ^  ^ 

and 


Ik. 


M 


N 


001-030 


CO, 


+  k. 


N. 


M  -  k  C02 

f 001-030  w2  ^001-030  “2 


CO.  ”  *f...  .?  \  +  kf 


He 
001-030 


N, 


He 


and  they  are  related  to  S?2  transitions. 


001-030 


is  the  rate  constant  that  a  CO^-molecule  in  quantum  level  (001) 
will  suffer  a  collision  with  another  C0o-molecule  and  make  a 
transition  to  level  (030)  per  second. 

Likewise  for  the  others. 
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Q 


and 


kf  ^  kf 


=  k. 


=  k 


■f  *s  related  to  SP.  transitions. 

020-021  1 


000-001  * 01 0-011 
I's  are  the  intensities  of  the  radiations  (refer  to  fig.  2) 

9  =  0(01  0)  *=  667.3cm-1:  The  energy  of  the  level  C02(0110) 

9  =  [0(001)  -  0(030)]  =  380m-1:  the  energy  difference  between  the  levels 
C02(001)  and  C02(030).  0 

We  note  that  In  our  computer  program  these  equations  have  been  divided 
through  Xk^M,  and  the  nondimensional  time,  t'  has  been  defined  as: 

t'  =  tZk*1  M 
fl-i-0 


o 


u 
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APPENDIX  C 


Four  Tempera Lare  Computer  Program 

Introduction: 

A  new  and  simplified  four  temperature  model  computer  program  was 
developed  to  explore  the  total  pressure  losses  of  various  geometries  and 
mixtures  and  also  to  support  gain  measurements  of  both  the  high  temperature 
studies  and  the  supersonic  nozzle  development. 

The  new  program  uses  an  approach  analogous  to  the  previous  two  temperature 
program  (1),  but  it  includes  improvements  other  than  using  four  temperatures 
instead  of  two.  As  before,  the  program  calculates  the  parameters  (p,  T,  gain, 
etc.)  along  any  given  streamtube.  In  addition,  the  new  program  includes 
Doppler  broadening  effects  on  the  line  center  gain  coefficient,  entropy 
production,  and  a  calculation  of  the  total  pressure  and  temperature. 

It  avoids  using  y  (=  cp/cv)  in  the  flow  equations  and,  with  four  temperatures 
has  a  more  realistic  lasing  routine.  The  effects  of  small  percentages  of  H20, 
c0»  °2,  and  Ar  can  also  be  included  in  this  new  program. 

Flow  Equations: 

A  set  of  flow  equations  has  been  derived  which  are  more  practical  for  this 

application.  A  necessary  result  of  considering  vibrational  energy  is  that  y, 

the  specific  heat  ratio,  varies.  To  avoid  using  y,  there  has  been  defined  a 

variable  B  to  replace  M,  the  Mach  number,  where  B2  =  —  =  yM2. 

RT 

The  quasi-one-dimensional-inviscid  flow  equations  were  developed  in  terms 

of  B  and  in  a  stepwise  formulation  for  computer  application.  Stepwise  area 
A2 

ratio,  is  the  only  "known"  of  the  equations.  The  momentum  equation  includes 
1  P  +  P 

the  wall  pressure  term  from  area  change,  (-^  ^(A^).  The  energy  equation 

‘S  for  translation  and  rotation  only,  since  Q  is  calculated  from  the  vibra¬ 
tional  energy  relaxation  equations. 


Continuity: 


(1) 


Zi  ti.LlS 

p.  a2mt, 


Momentum: 

Energy: 


TT  “  (1  +  7T  +  2  Bi)  !  (1  +  a1  +  2  a1  B2J 

1  :  i 

C  Q 

{  R^)T,R  +  h  Bj)  Tj  +  Q  =  {  -|)T>R  +  %  B2}  Tj 


(2) 


(3) 


where  g-)T>R  includes  only  the  translational  and  rotational  degrees  of  freedom, 
and  Q  is  the  heat  transferred  into  translation  and  rotation  from  vibration. 


Vibrational  Relaxation  Equations: 

The  visrational  relaxation  equations  were  written  f.or  four  temperatures: 

T,  translational-rotational;  l’L,  the  lower  modes  of  C02  (symmetric  stretch 
(1)  and  bending  (2));  Ty,  the  upper  mode  of  C02  (asymmetric  stretch);  and  T 
the  N2  vibrational  temperature  (see  Figure  (1)).  The  corresponding  energies  in 
modes  are  EL(=E^  +  2E2),  E^  and  E^.  The  assumption  of  equilibrium  within 
a  mode  is  made,  and  the  bending  mode  is  assumed  to  be  in  equilibrium  with  the 

symmetric  stretch  mode.  There  is  no  assumption  that  TL  is  "nearly  in  equili¬ 
brium"  with  T. 

In  order  to  use  the  technique  of  Stollery  and  Smith (2),  the  relaxation 
equations  were  derived  into  the  Bethe-Teller  form: 


dE 

dt 


J  (E*  -  E} 


where  E  is  the  vibrational  energy,  E*  is  the  equilibrium  vibrational  energy, 
and  T  the  relaxation  time  (  a  function  of  p,T).  In  order  to  do  this  without 
approximations,  some  sloyly  varying  quantities,  all  nearly  unity,  were  intro¬ 
duced  (ZA,  Zg,  Zc).  The  relaxation  equations  became: 


2 
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Figure  1.  Vibrational  Energy-Level  Diagram  for  C0„  -  N„  Syst 


Tiacre 


V  {E"  “  N  w  *  *  v 


Ej  E  1 

ZA  -  if*  Ef 

a3  /!_  _  1  v 

exP  fcr*  rf*  (0.4-0.)  '■T  1  T^)  -  1 

B  0j— 

<=*P  (j-)  -  1 
®3  1 E 

ZC  "  e~  ZD  +  T~  0  (ZD  ~  X)  » 

N  yC02  N 


®N  "  03 
;  Zp  =  exp  - - i) 


h  "  Ei  +  2Ez 


^CO  61 


E1 " 


exp  (— )  -  l 
L 


;  e2(tl)>  EpCTp).  ENaN)  are 'analogous 


^CO  i  =  mole  fractions  1 

i 

0l»  02»  03’  0N  ■"  characteristic  vibrational  temperatures 
note:  q  -  ^  +  y  +  y  -  a-  +  ^  +  e^)  '  ■ 

TI.e  first  term  In  (4)  represents  the  deectivatlon  of  the  lower  levels  to  the 

ground  state.  The  last  term  In  (4)'  and  the  first  In  (5)  represent  the  transfer 

of  energy  from  the  upper  level  to  the  lower  levels.  The  last  term  in  (5)  and 

the  first  of  (6)  represent  the  coupling  of  the  upper  level  and  the  N.,  vibration. 

The  last  tern  of  (7)  is  the  slow  deactivation  of  the  »2  vibration  directly  to 
the  ground  state. 

To  solve  the  relaxation  equation,,  the  equations  have  been  written  in 
difference  equation  form  suggested  by  Stollery  and  Smith<2).  Note  that  (5 
and  (6)  can  he  solved  independently  of  (4).  By  solving  for  E,  ts  ?  function  of 
V  ^  is  solved  without  iteration.  i  1 
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} 


StMMrataS’ttCi  si 


o 


lasing  is  presently  accomplished  by  extracting  all  available  energy  in 
the  J02  upper  level,  while  appropriately  raising  the  lower  level  energy  until 
zero  gain  is  reached.  This  corresponds  to  an  infinite  laser  intensity. 
Relaxation  Time  Equations: 

The  relaxation  times  V  V  V  and  (t^  are  from  Taylor  and 
B.tterman's(  data.  The  remaining  T^  are  frna’the  data  of  Rosser,  Wood,  and 

(4^  -n/o 


Gerry  .  The  pT  versus  T  1  3  equations  are  approximations  to  the  best  fit 

curves  given  and  are  in  agreement  with,  and  fit  the  data  better  than,  those  of 
Munjee^  and  Anderson^. 

i 

The  equation  used  to  find  the  relaxation  time  of  a  gas  mixture  from  its 
components  is: 


* 


pT  all  species  (pr) 


species 


species 


The  curve  fits  for  Ptl>  in  atm-sec: 


300  -  460eK  (PTL)C02  *  2 (10“6)  exp(24T“1/3) 

(pVn  *  9<10"7>  (24.7T"1/3) 

i  2 

(pTL)He  =  9. 34(10~8)  exp(8.96T"l/3) 

460  -  2000°K  (PTL)C02  =  3. 06 CIO"9)  exp(56.1T-l/3) 

r“l/3N 


(pVn2  “  1-6K10"8)  exp(55.8T~1/3) 
(pTL)He  =  1‘8  <10"8)  exp(21.6T-1/3) 
300  -  2000°K  (ptl)h  0  =  1.04(10~5)  exp(-55.5T~1/r3) 


The  curve  fits  for  px^,  in  atm-  sec: 


300  -  380°K  (PTu)CQ2  =  2. 04 (10” 6 )  exp(4.28T-1/3) 
380  -  2000°K  (pTu)  -  2.67 (10_1 °)  exp (68. 9T"1/3) 


300  -  615°K  (pVn2  =  6.  75 (10~8)  exp (34. 2T~1/3) 

615  -  2000°K  (pTy)^  =  4.3  (10"12)  exp (116. ST"1'3) 
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300  -  2000-K  <PVHc  -  3.5  (ptu)k 

(pTU}H  0  *  2i74  <10”7)  exp(-12.8T“l/3) 

2 

The  curve  fit  for  pT^,  in  atm-sec : 

300  -  1000°K  (ptN)N2  =  1.4  (10"11)!1*5 

1000  -  1400°K  (PTn)n  =4.5  (10~7) 

1400  -  2000°K  (ptN)N  =  3.5  (10~6)  expt--^) 

The  curve  fits  for  pt  ,  in  atm-sec: 

s 

300  -  2000°K  (pt  )  =  (ot  ) 

s'c02  ^pTs;N2 

<pTs)N  -  9.52  (10~17)  T  exp(278T“1/3) 

2 

<pVHe  =  2,88  <10'X3>  T  exp(lllT"1/3) 

<pTs>H20  *  7*6  (10-U)  exp(88T_1^3) 

For  small  fractions  of  CO,  02  or  Ar,  a  good  guess  is  to  use  (pr)  in  all 
cases. 

Summary: 

This  four  temperature  C02~N2  kinetics  program  is  quite  inexpensive,  yet 
retains  good  accuracy  (T  within  10°K  of  the  more  costly  programs).  It  also 
remains  flexible  and  can  handle  any  nozzle  flow  problem  with  vibrational 
relaxation.  The  program  is  being  expanded  to  handle  the  effect  of  varying 
laser  intensities.  Vibration-vibration  coupling  of  CO  with  C02  will  be 
included  when  the  data  becomes  available. 
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